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ABSTRACT 


Non Linear Internal Waves (NLIW) are ubiquitous and appear wherever a proper 
combination of stratified water, current and bathymetry occur. In recent years, they have 
also been proven of primary interest for acoustic oceanography, since they are known to 
play an important role in sound speed fluctuations in shallow waters. The predictability of 
acoustic variability caused by these waves has been somewhat limited, and largely based 
on direct numerical simulation. Through this project, we present a simple and 
computationally efficient analytic model based on coupled mode theory. For this purpose, 
narrow and broadband acoustic normal mode fluctuations in the 2001 South China Sea 
ASIAEX experiment are first examined. Then ASIAEX environmental data were used to 
characterize the space/time scales of NLIW sound speed structure. Einally, a comparison 
is conducted between the observed normal mode variability and the predictions from an 
analytic model utilizing the observed NLIW structure, such as its width, amplitude and 
speed. As the model proves to be very sensitive to these characteristics, an overview of 
the Synthetic Aperture Radar capabilities to retrieve them is also conducted. In this latter 
part of the project, particular attention is given to wind effects on SAR data. 
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I. 


INTRODUCTION 


With the emergence of asymmetric threats, modern naval warfare has moved from 
the deep ocean to the shallow waters, a fact that has necessitated a thorough 
understanding of this complex and variable environment. With regards to acoustics, much 
attention has been paid to the variability of the bottom geology and geoacoustic 
parameters due to the importance of bottom interaction in shallow water acoustics. 
However, it has also been recognized in recent years that the water column and its 
fluctuations also have a large impact on modulating acoustic signal in shallow water 
(Headrick et ah, 2000, Duda and Preisig, 1999). One aspect of this variability includes 
non linear internal waves (NLIWs), the so-called solitons or Internal Solitary Waves 
(ISW). These waves are accompanied by huge vertical displacements of water, upward or 
downward in accordance with the local density profile, and they affect significantly 
sound speed profiles. Thus, the present thesis aims at improving the understanding and 
prediction of the influence of solitons on sound propagation in shallow waters. Toward 
this end, this thesis tests a Dyson series-based method presented by Colosi (2007) using 
in-situ data collected during the Asian Seas International Acoustic Experiment 
(ASIAEX), conducted in May 2001 in the Northeastern South China Sea. Because of the 
importance of the inner characteristics of solitons in this theory, such as their half-width 
or their amplitude, this paper will also present an overview of the remote sensing 
techniques used to detect and characterize them. Particular attention was given to the 
wind effects while retrieving this information. 

A. OVERVIEW ON NON LINEAR INTERNAL WAVES 

Apel et al. (2007) have provided a thorough review of non linear internal waves 
and solitons, and this section presents only some of the main characteristics. 

As their name implies, internal waves travel within the interior of the ocean and 
are particularly characterized by strong vertical anomalies in temperature. As illustrated 
in Eigure 1, non linear internal waves are ubiquitous, appearing wherever the proper 
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combination of stratified waters, currents, and bathymetry occur. This characteristic 
makes them of primary interest from a military perspective. 


Solitons tend to oeeur in paekets, usually rank-ordered, with the largest 
oscillations appearing at the paeket front. They generally preeede an oseillating 
undulatory bore assoeiated with a depression of the equilibrium pyenoeline. Despite the 
growing interest on this topie, their generation is still ill monitored; however, hydraulic 
control by seamounts or continental shelves during energetic tidal cycles seems to be an 
adequate explanation in many cases (Ramp, 2006, personal communication). However, 
their physieal eharacteristic show significant modulations on diurnal, fortnightly, 
seasonal and semiannual time scales, and these changes are greatly associated with the 
ehanges in the water column stratification. Furthermore, non linear internal waves are 
evolving events. Aetually, as they propagate, their shape shifts from a single, relatively 
sharp peak to a series of spikes, one oseillation being added per Brunt-Vaisala frequeney. 
They also tend to slow down before finally being dissipated by different meehanisms, 
among whieh are radial spreading, bottom interaetion, instability and fluid turbulenee 
(Apel et ah, 2002). The lifetime of a soliton is on the order of a few days in the open sea 
to a day on the eontinental shelf. 
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Figure 1. Locations of internal wave imagery and data presented in Internal solitons 
in the ocean and their effect on underwater sound (Apel et ah, 2007). 
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The first scientists to give mathematical representation to this phenomenon were 
Korteweg and De Vries (KdV) in 1897 who based their work on the observations made 
by Russel in 1838 and 1844. Since this original work, many additional solitary wave 
equations have been derived, and new solutions, such as the dnoidal representations, have 
been found. However, the KdV equation or the extended KdV equation (Apel et ah, 
2007) offer results whose accuracy is sufficient for the purpose of this project and will be 
used to characterize non linear internal waves. 

Solitons can be described by a product of solutions to the weakly non linear KdV 
equation, D{x,t), and the eigenfunctions of the Taylor-Goldstein equation for the 

vertical structure function W); „(z): 


*‘max 

z,t) = Ya .IT, „ (z).D„ (x, t) 


( 1 . 1 ) 


D{x,t), the normalized vertical displacement of an isopycnal surface from its 
equilibrium level, andlT, „(z) are the solutions of the two following equations; 


dP 

dt 


- + c„ 


dP ^ dP d^P 
— + a.P .— + y —^ 
dx dx dx^ 


= 0 


( 1 . 2 ) 


d^W, 


k,n . 1 2 


dz" 


+ k\ 


N\z)-o/ 


1 rl 

a -f 


■W,Az) = 0 


(1.3) 


In the three previous equations,«, / and Cg are environmental coefficients 

describing the nonlinearity, the dispersion and the phase speed associated with the waves, 
and are functions of the reduced gravity and the thicknesses of the upper and lower layers 
as illustrated with the following equations; / represents the Coriolis parameter. Thus, if 

/ij is the thickness of the upper layer and ^2 the thickness of the lower layer, then a , jd 
and Cg are defined by: 

^ (1.4) 

2 /ij ./ij /?2 -hi + /?[ ./I 2 
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(1.5) 


p— ^0 A-^2 P\ a Pi -^2 

6 'p^\+p^.h^ 

^ _^ gKKiP2-Pi) ^ 1/2 

Pi A A 'A 


( 1 . 6 ) 


A^(z)is the buoyancy frequency, k characterizes the horizontal wavenumber of 
the soliton and (uis the radial frequency. Two boundary conditions finally add to make 
this set of equations a Sturm-Liouville system. The description of the physics behind the 
ISWs using a two-layer model is convenient since it offers a simple interpretation of the 
problem, nevertheless, any of the three previous parameters a,/ andcg, as well as the 

ISW mode functions can currently be calculated using a continuous stratification of the 
water column. In this thesis, both approaches are used, a continuous stratification for the 
section dedicated to the acoustics fluctuations associated with the solitons and the two- 
layer model during manipulation of SAR images. 

If the shallow water approximation holds and if the upper layer depth is small 
compared to the width of the soliton, the typical solution for a single soliton pulse is then 
given by: 

;7(v,0 = ?7oSec/zA^-^) (1-7) 

A 


where A = 



is the characteristic width and V =Cg 



the so called non linear 


speed. Figure 2 represents r/(x,t) as described by equation 1.7 whereas Figure 3 shows 
the dnoidal model, a more realistic but more complex solution to equations 1.2 and 1.3. 
ISWs are often seen as depression associated with downward vertical displacements of 
surface waters, elevation or upward water displacements are also observable particularly 
in very shallow waters where /Zj can be greater than /Zj. 


In this report, equation 1.7 operates for//(v,?) in the section dedicated to the 
remote parameterization of solitons using SAR systems, but is approximated by a 
Gaussian profde during work with solitons and acoustics in shallow water (Young, 2006). 
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Figure 2. Classic secant profile for KdV soliton at two different times, viewed in a 
moving coordinate system. Adapted from Internal solitons in the ocean and their 
ejfect on underwater sound (Apel et al, 2007). 
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Figure 3. Dnoidal model evolution over one diurnal period. This model includes the 
vertical structure function and attenuation, and adds one cycle every buoyancy 
period. Adapted from Internal solitons in the ocean and their effect on underwater 

sound (Apel et ah, 2007). 


5 




B. THE ASIAEX PROGRAM 


To support all the results presented here, this project used data collected in the 
South China Sea (SCS) during the second leg of the Asian Seas International Acoustics 
Experiment (ASIAEX) which took place in April and May 2001. The goal of the SCS 
project was to understand acoustic propagation through shallow water when strong 
variability in the form of fronts, eddies, and internal waves, is present (Ramp et ah, 
2004). In essence, the project was a progression from previous similar coupled physical 
oceanography and acoustics experiments, namely the Shallow Water Acoustic Random 
Media (SWARM) (Headrick et al. (2000)) and Shelfbreak PRIMER (Eynch et al. (1997)) 
experiments. As an added benefit, this experiment presented a unique opportunity for 
international cooperation in oceanography and a significant effort to improve the 
knowledge of non linear internal waves that occur significantly in this strategic area. 
Numerous works were published about this experiment, particularly in volume 29, 
number 4 of the IEEE Journal of Oceanic Engineering published in October 2004. The 
current section presents only the data used for this project. 



* S8 


Eigure 4. Map showing the locations of the ASIAEX moorings. The oceanographic 
moorings are indicated by stars labeled S2-S8. The shaded red boxes indicate 
acoustic assets, with an “S” indicating sources and “E-Array” the receiving 
arrays. Adapted from Internal Solitons in the Northeastern South China Sea. Part 
I: Sources and Deep Water Propagation (Ramp et ah, 2004). 
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In order to get a elose look at the effeets of solitons on sound propagation in this 
shallow and variable environment, measurements of water eolumn properties but also 
vertieal and horizontal distributions of aeoustie pressure were eolleeted using numerous 
moorings (Figure 4). To resolve the high variability of the medium, temperature, salinity 
and pressure were eolleeted at high resolution in spaee and time. In partieular, the array 
eonsisted of seven densely instrumented oeeanographie moorings deployed aeross the 
eontinental shelf. With this eonfiguration, the moorings spanned depths from 800m to 
72m, along the primary aeoustie paths, reeording at a 2 minute-sampling rate (Ramp et al, 
2004). During this experiment. Ramp et ah, deseribed earefully the generation and 
propagation of the solitons under observation and partieularly noted the evolution of the 
solitons as they moved shoreward. Aetually, Ramp et ah, diseemed that the extended 
KdV model gave better results than the basie KdV theory in a deseription of the solitons 
as they shoal. The former theory eannot adequately aeeount for the ehanges in wave 
shape and speed. Numerous results eoneerning the amplitude, speed and width were also 
eolleeted. 

Regarding sound transmission, one sound souree operating at 224 and 400Hz was 
deployed on the 350m-isobath. The souree transmitted phase-modulated signals at a 
earrier frequeney of 400 Hz with a bandwidth of 100 Hz and a souree level of 180 dB (1 
pPa). Frequeneies on the order of a few hundreds hertz are aetually known to limit 
seattering and bottom effects. The phase modulation employed was a 5.1 l-s-long511- 
digit pseudorandom sequence resulting in a compressed pulse of 10 ms-resolution after 
matched filtering. The acoustic data were received by an L-shaped hydrophone array, 16 
were moored vertically in the water column and 32 hydrophones horizontally. The 
sampling was realized at a rate of 3.2 kHz over the first three weeks of May 2001 (Chiu 
et ah, 2004). This work focused only on data recorded by the vertical array (VLA) for the 
400Hz source. 

Many significant results arose from these different observations. For instance, 
Chiu and al. (2004) highlighted the fact that the variability of the sound speed occurred 
within three distinct spectral domains which they labeled tidal, subtidal or mesoscale, and 
supertidal, that is due to internal waves. They also discovered that this variability could 
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disturb the acoustic field by a few decibels due to the vertical redistribution of acoustic 
energy, particularly in the presence of internal waves. This phenomenon is known as 
mode coupling; it was already observed during the SWARM experiment and is within the 
purviews of this work as well. 

Concomitantly with this experiment, SAR images of the area were taken and will 
allow us to look at the relevance of this system in detecting and characterizing solitons. 
During the time frame of the experiment, only images from the Canadian system 
RADARS AT were available; these will also be used in this project. 

C. STATEMENT OF RESEARCH 

1. Acoustics 

The present analysis aims to study the acoustic intensity fluctuations, such as time 
modulation and spreading, caused by the presence of solitons along the acoustic path. For 
this work, acoustic pressure and energy are described in terms of normal modes. This 
decomposition originates from the resolution of the Helmholtz equation using a 
separation of variables technique. The natural frequencies governing this decomposition 
define mode functions and are highly dependent on the boundary conditions. This 
common analysis of the propagation of acoustic signals was initially designed for range 
independent environment but can also be generalized to range dependent medium. 
However, in these latter kinds of environment, normal modes are generally no longer 
independent, which may conduct to the vertical redistribution of acoustic energy. Three 
main points of views are commonly expressed regarding this topic: adiabatic 
approximation, weak mode coupling and strong mode coupling. With the first 
approximation, interaction between modes is neglected, this assumption was proven to be 
erroneous regarding this topic (Rousseff et ah, 2002). The second approximation assumes 
that the primary interaction occurs between consecutive mode orders and is used all along 
this thesis. Finally, the third approximation is generally associated with a complex 
bathymetry where the bottom effects exceed the effects caused by the water column. In 
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this work, we also looked at this redistribution of aeoustie energy by analyzing the 
exehanges whieh oeeurred between modes during the Asian experiment for different 
environmental conditions. 

The second goal of this thesis was to test the ability of a simple and fast analytic 
model to predict the changes induced by the ISWs. This model designed by Colosi (2007) 
builds upon works from Dozier and Tapper! (1978), Dozier (1982), and Creamer (1996). 
The model also clarifies the work of Zhou et ah, (1991), which noted the resonance 
condition that exists between various wavenumbers of a NLIW and the beat 
wavenumbers of the acoustic normal modes. The basic solutions of this problem can take 
the form of the Dyson series from quantum mechanical perturbation theory (Sakurai 
(1985)). Indeed, they are designed to solve for the change in modal amplitudes after 
acoustic propagation through ocean NLIWs. This treatment is confined to coupling in the 
depth-range space, and thus azimuthal coupling is ignored. Moreover, as this theory only 
assumes weak single scattering, it is limited to frequencies between 100 and 400Hz and 
breaks down at higher frequencies. Finally, this approach can treat either single NLIW or 
packets of NLIWs. 

Young (2006) verified the relevance of this theory by successfully testing it 
against numerical simulation to determine if the total acoustic pressure field could be 
accurately predicted in an ocean environment with ISWs. The present study aims to 
validate the theory in a real environment, based on the oceanographic and acoustic data 
provided by the ASIAEX experiment. 

2. Remote Sensing using SAR Images 

Synthetic Aperture Radars (SAR) are radar systems which have improved their 
resolution by extracting information from the spectral shifts caused by Doppler effects. 
Because of their global coverage and their specifications, SAR satellites act as primary 
providers for information about NLIWs. As it will be detailed in the fourth chapter of the 
present document, SAR satellites possess a unique ability to detect these waves by using 
the interaction between wind waves and the current induced by ISWs (Alpers (1985)). 
They also allow users to retrieve some essential soliton features such as their half-width, 
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their wavelength or their amplitude (Liu et al. (2004), Zheng et al. (2001)). These latter 
eharaeteristies are essential for the previous aeoustie eonsiderations. 

The goal of the fourth and fifth seetions is to synthesize the main results of the 
works made in the last 20 years regarding soliton deteetion by SAR systems. Thus, these 
seetions will introduee basie theories about soliton deteetion by SAR satellites. Attention 
will partieularly foeus on two previous works. The first, a study by Zheng et al. (2001), 
presents an interesting method of retrieving soliton information from SAR images. The 
seeond, a study by Liu et al. (2004), gathers the results observed using SAR images 
during the ASIAEX experiment. As mentioned previously, this work also earefully 
investigated the wind effeets to determine how they ean bias the soliton information. 
Unfortunately, beeause of the laek of data, the study will offer only a eouple of images. 
Thus, if the statistieal aspeet of the study does not prove very eonelusive, the projeet still 
provides an interesting qualitative overview of the advantages and limitations of SAR 
satellites with regard to this topie. 
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II. THEORY AND PRESENTATION OE THE PROCESS 


A. THEORETICAL MODEL 

1. Shallow Water Acoustics 

The present model was designed to compute the acoustic intensity fluctuations 
induced by NLIWs in shallow waters. The model assumes only single weak scattering, 
and thus provides noteworthy results solely for frequencies ranging from 100 Hz to 
400Hz, which are of primary interest for submariners. 

Regarding the spectral domain and the distances used, the attenuation of the 
signal due to the water column is negligible; only the attenuation due to the bottom is 
considered. Moreover, the model considers that the total acoustic wavenumber will have 
a small imaginary part: 

k{r,z,co)= ^ +i.a{z). (2.1) 

c{r,z) 

Here a{z) is the attenuation coefficient and is considered non-zero in the seabed. 
c(r, z) is the sound speed profile at a distance r from the source and at depth z . k and co 
are respectively the wavenumber and the radial frequency of the signal. 

The present theory furthers assumes that the incoming soliton will convey an 
isolated sound speed perturbation to the background sound speed profile. In their study, 
Chiu et al. (2004) noticed the importance of the tidal and subtidal variabilities; however, 
the present study sets these effects aside to focus only on the high frequency variability 
induced by non linear internal waves. With this assumption, one can decompose the 
sound speed profile, such as: 

c(r,z) = c(z) + Sc(r,z) (2.2) 

where c(z)is the background sound speed profile and Sc{r,z)is a perturbation due to 
some NLIWs. 
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The development outlined by Dozier and Tappert(1978), Dozier (1982), and 
Creamer (1996) deeomposes the acoustie pressure field for the radial frequeney coin 
terms of the aeoustie normal modes of the range independent problem sueh that: 

N A (f.\ 

P(r,z) = ^^)=^.(p„(z)- (2.3) 

n=i yjr 


A^(r) represents the amplitude and eharaeterizes the range-dependent effeets of 
seattering on aeoustie pressure. A^{r) is a solution of equation 2.6 with A„(0)being the 
mode amplitude at the souree. ^z!'^(z)is the n* normal mode funetion and is a solution of 
equation 2.4. Aetually, ^z)^(z) represent the unperturbed normal modes ealeulated in the 
absenee of ISW and all the variability eonveyed by the ISWs is ineluded in the amplitude 
term A„(r): 


[cr(z) 


5 1 


5 . 172 


dz a{z) dz 




(2.4) 


where 


k{z) = 


—— + i.a(z) and a(z) is the baekground density profile. The 
c(z) 


orthonormality relation assoeiated with ^„(z) is given by: 


° cr(z) 


<lz = 


In equation 2.4, , the eomplex wavenumber is written +/.«„, where is 

the usual modal wavenumber in the absenee of attenuation , and is the eomplex part 
of the wavenumber, obtained using the perturbation theory: 




61 r «(z).^„"(z) 
kn 0 c{z).a{z) 


dz « k„ 


(2.5) 


To eomplete this elose Sturm-Liouville problem, the boundary eondition at the 
surfaee is defined by ^zi„(z = 0) = 0. Then, in order to take into aeeount the bottom 
attenuation and to keep this model eomputationally effieient, the bottom boundary 


eondition is ehosen to be 




= 0 at some seabed depth far below the water-bottom 


z=h 
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interface, generally 3 or 4 times the bottom depth. Finally, to overcome computational 
difficulties related to the interface between the water column and the seabed, the 
discontinuous density or sound speed profiles were substituted by sharp but continuous 
profiles. 

In the quasi-static, narrow angle, weak forward scattering approximation Dozier 
and Tappert (1978), and Dozier (1982) show that the mode amplitudes can be 
approximated by the one-way wave equation; 



Mi=i 


( 2 . 6 ) 


where the symmetric coupling matrix , initiated by the perturbed sound speed 
profile, is given by: 


Pmnir)= 

iKK) 


00 

■\p{r,z). 




dzwith p{r,z) = 


dc{r,z) 


Here, kg = — and Cg is a reference sound speed for the local water column. In 
Co 

equation 2.7, the coupling matrix does not depend on the bathymetry and is only a 
function of the sound speed perturbation induced by the ISWs. It is one of the strongest 
assumptions of the model which assumes that the adiabatic approximation holds once the 
energy is redistributed among the modes by the solitons. 

Finally, it is useful to remove the rapid oscillations, which also makes the 
equation symmetrical, by introducing (p^ such as: 

(Pnir) = A„(r).exp(-/./„.r).k^'". (2.8) 

Then, equation (2.6) becomes: 

^ = -6^ Pn,n -c' '"" ■>„ with (2.9) 

or 
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2 . 


Solution using the Dyson Series 


Although it can be significant, the sound speed perturbation induced by a solhon 
is small compared to the background sound speed and generally less than 30 m/s. From 
this observation, it can be deduced that and p are small, as expected from equation 
2.7, and that the mode coupling is relatively weak. In fact, this strong assumption 
suggests that the problem can be solved using a perturbation theory, such as may be 
found in standard quantum mechanics texts. Actually, the Dyson series proved valuable 
for this purpose and were used to solve this set of equations. Through this approach, it 
can be shown that: 

ft (S) = ft (0) - i.2; (0).J {r'U'-'-dr’ (2.10) 

m=l 0 

with ^„(0) = kyi|p(0,z).-^^Jz. (2.11) 

^„(0), like A^(0) earlier, represents the amplitude of the normal modes at the 
source and is greatly dependent on the depth of the source. Using the Fourier transform of 
, equation 2.10 becomes: 

TV 00 

(Pni^) = J krnik )]^ Jr’.Jk (2.12) 

m=l ^ 0 


1 ^ 

with p^{k) = — • J Pmnirle' '-'dr. 


Moreover, for the considered range of frequency, as R is large and a is small, 

R 

one can approximate //„,„(/?,k) = Jr' by a Dirac function such as 

0 

H„,„(R,k) ~ 27r.S(k^„ -k) where k „ =k^-k„ is the beat wavenumber between mode m 
and n. The Dirac representation of H^^{R,k) is really a cornerstone of this model since it 

enhances the fact that the resonant coupling condition is conditioned by the match 

between the beat wavenumber and the horizontal wavenumber of the ISW. 
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Finally, the equation beeomes: 


(Pn (^) = (0) - i-X (^) 

m=l 


with7^„(i?)= J p^^{k).H^^{R,k)dk~27t.p„„{Kn) 


(2.13) 


and, A.(^) = [A.(0)-i.|;(^)‘'lA,(0)J,.(/f)].<."-'. (2.14) 

m=l % 

In this equation, J is the seattering matrix for mode m and n. For multiple 
solitons, beeause of the linear relationship between ^^(i?)and p^„{k), will simply 

M, 00 N, 

beeomes. where is the number 

-00 7 =1 

of solitons. 


3. Non Linear Internal Wave Sound Speed Perturbations 


In the previous set of equations, the sound speed perturbation ereated by the 

Sc(r z) 

presenee of one or more solitons is ineluded in the variable p{r, z) =-. Munk and 


Zaehariasen (1976) derived the following equation for p ineluding the internal wave 
vertieal displaeement; 


= (2.15) 


Cq dz 


Cq ' dz 


In this equation, (—) represents the potential sound speed gradient whieh is used 
dz 

beeause internal waves move water pareels adiabatieally. The shape of the vertieal 
displaeement indueed by a soliton was already presented in seetion LA and is redueed to 
the first mode of equation 1.1 as it is generally assumed in the literature. F(r) is 
essentially deseribed by equation 1.7 but is approximated with a Gaussian soliton in this 
seetion as deseribed by Young (2006). It is important to note that ehanges in the NLIW as 
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it propagates, specifically its bore-like structure, are not considered; these changes 
introduce analytical complexities which encumber the computational efficiency of the 
model. However, the linear decay of the waves with the bathymetry is simulated. Finally, 
one can rewrite equation 2.7 thus: 


Pmnir) = 


iKK) 


^]p{r,z)' 




m n / 0 


cr{z) 


dz = C^„.F{r) 


withC “ 


Co iKK) 0 dz Cj{z) 


(2.16) 

(2.17) 


Here matrix characterizes the contribution to coupling due to the vertical 
structure of the acoustic modes and of the soliton. As such, C is determined by the 
overlap between the depth structure of the NLIW and modes n and m. For a packet of 
Gaussian solitons located at the distance rQ{j) from the source, one can rewrite such 
that: 


on 


= with 




-+i-r„(]).t 


(2.18) 


y=i 

As explained in section II.A.2, can be viewed as a Dirac function centered 
, and then, using equations 2.15 and 2.18, the equation becomes: 


4W = [A,(0)->'.'SZi(^)‘"A,(0)c,„(i).A(;> 

y=l "1=1 


■A{A 


+‘-ro(j).l 


].e"" * . (2.19) 


The latter equation quantifies the mode energy for the n**' mode received at a 
distance R from the source when a group of ISWs are located at the distances r^ij). This 
equation can be considered as time dependent, since it oscillates with the displacement of 
the ISW whose position is given by rg( 7 ), and is then a function of time. Moreover, as 

(R) is still a complex number, it contains both amplitude and phase information. Later, 
this analysis will involve separation of these two quantities and introduction of the 
expressions for the mode energy and the time spread. 
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Equation 2.19 also emphasizes the role of the main parameters involved in this 
problem, partieularly A^(0), the amplitude at the source which sets an average value, 

A(j)the width of the j* ISW and the beat wavenumber which act together to 
modulate the oscillations. Actually, plays a significant role for neighbor modes only. 
Finally the coupling matrix C^^{j) clearly appears as a centerpiece of the perturbation. 

B. ANALYSIS 

1. Preliminaries 

As introduced earlier, the goal of this project consisted in applying this theory to 
in-situ data in order to check its relevance for operational purposes. Hence, all the 
acoustic and oceanographic data involved in the numerous computations presented in this 
work were provided by the acoustic array described in section 1. The local temperature 
was measured at eight different levels at S4 and S5 and at 16 levels at S7, while salinity 
measurements were respectively collected at three, five and four different levels at the 
three locations. Along the vertical array, 16 hydrophones expanded from the 42m-depth 
to the bottom, located 120m below the surface, and recorded acoustic pressure at 3.2 kHz. 

Based on the temperature variations in time and depth at these three locations, 
three periods of time were thought to be of first interest for this analysis. These simple 
temporal windows were ultimately chosen because they were highly characteristic of a 
fluctuating environment. Thus, a first sequence with a single soliton, a second one with a 
well-defined soliton packet and finally, a third one presenting no soliton activity were 
used. The last sequence is taken as a reference for periods with very little acoustic 
fluctuation. The first timeframe extends from May 5, 2001 around 1945GMT to May 6, 
2001 around OlOOGMT. During that period, a single soliton propagated shoreward along 
the acoustic path, its characteristic half-width evolved from 240m near S7 to 290m near 
S4, its speed was estimated to be around 2 m/s between S7 and S5 and slow down to 1 
m/s between S5 and S4, its amplitude also varied from 65 to 40 m. Compared to the 
results presented by Duda et al. (2004) or Ramp et al. (2004), this soliton seems to be 
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much faster between S7 and S5 but this latter result is only an artifaet. The direction in 
which the soliton propagated was not parallel to the acoustic path; therefore, the pulse did 
not cross perpendicularly with the soliton. Then, with respeet to the pulse, the wave 
appeared larger and seemed to propagate faster. Moreover, the decrease in speed after S5 
was probably due to the interaetion between the soliton and the bottom. It was 
accompanied by a horizontal spreading and a vertical squeezing of the soliton which 
refracted as it shoaled. 

The second sequence occurred on May 6, 2001 between 0700GMT and 
1500GMT. The soliton packet was composed of three solitons, the same 
spreading/squeezing phenomenon and decrease in speed were observed as the solitons 
moved to the coast. In that case, the initial amplitude of the leading soliton was 110m 
and deereased to 60m as it passed S4. The eharaeteristic half-widths varied from 320m to 
450m. The third environment was observed on May 7 between 0330GMT and 0730GMT. 
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Figure 5. Temperature field from mooring S7 (x=0) to mooring S4 (x=33km) on 

May 5, 2001 at2306GMT. 

It is noted here that the observed values for characteristic width or speed at S7 and 
S5 still agree quite well with the theory described in section I which is not true anymore 

when the wave interacts with the bottom between S5 and S4. 
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Figure 6. Temperature field from mooring S7 (x=0) to mooring S4 (x=33km) on 

May 6, 2001 at 0920GMT. 
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Figure 7. Temperature field from mooring S7 (x=0) to mooring S4 (x=33km) on 

May 7, 2001 at 0650GMT. 
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Once these three periods of time were identified, the seeond step of the project 
consisted in determining the assoeiated baekground sound speed profiles which would 
serve as input in the theoretical model. These baekground profiles were based on the 
recorded temperature, pressure and salinity measured at S7 and S4. In order to remove 
the effect of the non linear internal waves in these ealculations, a low pass filter was also 
applied to the data. As illustrated in Figure 8, the profile did not change significantly 
from one time to another, or between S7 and S4. A slight deerease in sound speed, 
around 2m/s, was observed between S7 and S4 in the presence of solitons. This effect 
may have been caused by the bore-like trail of the soliton on the measurements at 
mooring S7, which was hit first by the waves. 
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Figure 8. Baekground sound speed profiles for the three periods of time; quiet time, 
1®‘ soliton time (1) and 2"^* soliton time (paeket) at S7 (blue) and S4 (green). 
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The third step of this work eonsisted in ealeulating the normal modes assoeiated 
with these environments for a 400Hz-souroe with a lOOHz bandwidth. However, only the 
normal modes evaluated at S4 for a 400Hz-souroe were used in this theory. This 
ealeulation does not present a great eomputational ehallenge but does raise two important 
questions regarding the influenee of the bathymetry and the frequeney dependenee of the 
normal modes. The observations and assumptions resulting from these interrogations are 
presented in the next seetions. 

The final step of this projeet was the basie eomparison of the modal amplitudes 
ealeulated through the simple analytie model with those retrieved from the modal 
deeomposition of the aeoustie pressure, measured on the hydrophones as the soliton 
propagated between the souree and the VLA. All these results are presented in the third 
ehapter of this report. As the aeoustie pressure was reeorded on a diserete array of only 
16 deep hydrophones, the next seetion addresses the relevanee of this eomparison by 
analyzing the eross-talk matrix. 

2. Assumptions 

Young and Colosi (2006) tested the present theory against the simulation to 
determine if the total aeoustie pressure field ean be predieted in a flat bottom ISW 
environment. For this purpose, they used different frequeneies, souree depths, a 
theoretieal point souree and Gaussian solitons. Although their results were eonelusive, 
this present study will test the theory in a realistie environment with in-situ data and a 
variable bathymetry. With regard to simulation, eomparisons of theoretieal results to 
observed measurements present other ehallenges. Speeifieally, the diserete and relatively 
small number of hydrophones moored during ASIAEX 2001 raises a eentral question; is 
it possible to eompare objeetively simulated data from this theory with observed data 
eolleeted through only 16 hydrophones? Then, if the model works well for a point souree, 
one must also antieipate its behavior with regard to a pulse eovering a larger speetral 
band. 


21 



a. Frequency Dependence Mode Processing of ASIAEX Data 

From the normal mode theory as seen previously, for a source at a depth , the 
acoustic pressure field at the radial frequency nig and at a range R can be written such 
that: 


FfeO = (2-20) 

n=\ 

with . 


For this mode processing test in an unperturbed environment, this work considers 
only acoustic losses caused by attenuation due to the bottom and omit, for computational 
efficiency, cylindrical spreading and attenuation due to the water column. 


By applying the standard orthonormality relation, one then obtains: 




ik„(a>a)R _ 


piz) 


dz. 


( 2 . 21 ) 


For a pulse, the pressure will be given by: 

^ cOq+Xco 

p{z,t) = ^ j a„{o)).^„{o),z).e““do). ( 2 . 22 ) 

^^=1 oJq-Ao} 


The theory states that frequency plays a role in the definition of the normal modes 
at a particular location. However, as shown below in Figure 9, it appears that its influence 
is negligible compared to the amplitude of the modes for a lOOHz-spectral window. This 
statement can also be verified for any other modes and is generally truer near the bottom 
than near the surface. Thus, one can assume that (l)„{co,z) ~ (l)^{a)p,,z) at any depth and for 
any frequency. 
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Mode 1 



Figure 9. Comparison of first normal modes obtained at the loeation of the VLA 
(depth =125m) for different frequencies ranging from 350 FIz to 450 Hz. 


From the previous definitions, it is also apparent that the frequency influences the 
phase of the acoustic pressure field through the phase terme'*^"*^®^* . In order to account 

for the different phase speeds of the various frequencies and dispose of the phase term, 
one can introduce the modal group speed and make the following approximation; 

a\)+Aa) tyo+A<y ty^+Ao; 

j a„{o}).e-'“do}= j j (2.23) 

(Oq-ISlCO co^-ts.(o tyg—A<y 

coq+Aco coq+Aco ^ 

f f e-‘“^‘-‘"^dco where t„=R/—. (2.24) 

^ 3Jc 

coq-Aco cOq-Aco 
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In order to evaluate this result, one compares simulated pressure fields with and 
without the approximation. Figure 10 illustrates the relevance of one such calculation. 


Comparison between data 1 (no approximation) and data 2 (approximation) 



Figure 10. Comparison between the acoustic pressures calculated at a distance 
R=30km from the source (/o = 400//z - A/ = 20//z ) with (datal) and without 
(data2) any approximation of the amplitude and the phase term regarding 

frequency. 


Using (2.21), (2.23) and (2.24), the equation becomes; 






1 f 

p{z) 


(Oq+Is.co 

where (t) = j ^ 


(2.25) 


From these results, one can assume the possibility of amplitude retrieval for any 
modes based on the acoustic pressure field and the normal mode theory. Thus, evaluation 
of the model’s accuracy is easily achievable by comparing these values to those 
simulated by the model. 
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b. Cross-Talk Matrix 

Because of the discrete number of hydrophones, the orthonormality matrix, also 
called the cross talk matrix, is not the perfect identity and tends to skew results. For a 
discrete array, the orthonormality matrix is also approximated by; 


r 


mn 




(2.26) 


As it is represented in the following figure, the orthormality matrix is no longer 
perfect; the diagonal terms are not all equal to one and some of the non-diagonal terms 
are not null. Flowever, because only the first modes are of interest, this matrix still proves 
promising. Indeed, for the 5 first modes, the diagonal terms are still greater than 0.9 
which appears to be strong enough for the purpose of this study. Note that this theory 
would partly break for a shallow source since the energy would be concentrated in higher 
modes in that case. 




Figure 11. Cross-talk matrix for a sparse array of 16 hydrophones - the frequency is 
400FIz and the number of retained modes is 31 modes (trapped modes). 
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The last step of this preparatory work is to cheek the ability to retrieve the 
amplitude of the acoustic pressure at a distance R equal to 30 kilometers from the perfect 
source in an undisturbed environment and with a sparse array, through equations 2.27, 
2.28 and 2.29. 

(2-27) 


(^m(®o))me[l,6] ~ 

7=1 




ut) 


)mG[l,6] 


(2.28) 


16 




7=1 tv 


)me[l,6] (^m)me[l,6]-(rm«)me[l,6] ' (2.29) 

«g[1,6] 


Note; Equations 2.28 and 2.29 calculate the same coefficients but a perfect cross-talk 
matrix is assumed in equation 2.28 while the observed one is used in equation 2.29. 


Amplitude of the 6 first modes using different methods at t=19.74-44 s 



Figure 12. Comparison of the mode amplitude (amplitude only) recorded 30 km away 
from the source, and retrieved using three different methods. 
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Based on the results displayed in the above pieture, one may eonelude that 
equation (2.28) should be aeeurate enough to obtain the amplitude of the first six modes. 
From that, one may evaluate the relevanee of this theory by examining mode eoupling as 
the solitons appear. 
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III. RESULTS AND OBSERVATIONS (I) 


A. GENERAL COMMENTS ABOUT THE ENERGY RECEIVED AT THE 

VLA 

Using the acoustic pressure measurements recorded at the receiver, this study was 
able to retrieve the mode amplitudes for the first three modes, which were used to 
conduct the comparison. Indeed, this model assumes both weak coupling between modes 
and adiabatic propagation of the mode energy. Combined with the deep location of the 
source, these two assumptions make these modes the most energetic and relevant for this 
project. 

As they will appear frequently in the following figures and comments, it is also 
important to introduce the two time scales that frame the observations. The first one, t, 
corresponds to the geotime and sets the temporal reference of the experiment with regard 
to the Greenwich Mean Time (GMT). This time scale is somewhat discontinuous, since it 
only accounts for the beginning of each sequence and each pulse. In fact, the acoustic 
data observed at the array were not recorded continuously; rather sequences of 88 pulses 
were sent every 30 minutes from the source, 5.11 seconds separating the launch of two 
successive pulses. Then, only seven minutes of acoustic information were recorded per 
sequence and are available for this comparison. Any referenced time in Figure 13, 15 or 
17, such as 1815GMT or 1845GMT, is also a composition of 88 pulses, which are 
vertically displayed. Due to technical issues, it was not possible to retrieve the acoustic 
pressure for all the pulses in sequences such as 2345GMT or 0045 GMT in Figure 13; 
these sequences also appear narrower on the display. A second time scale, named r, was 
also introduced to analyze the data. This second temporal scale represents pulse travel 
time between the source and the receiver. 

Finally, this study chose to work with the mode acoustic energy, defined by: 

/„(r,r) = U.(r,r)P. (3.1) 
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Thus, Figures 13, 15 and 17 display the relative mode aeoustic energy for mode 1 
with respect to t and r for any of the three case studies presented in the second chapter. 
The energy is expressed in decibels in these three figures and represents the fluctuations 
in mode energy relative to the maximum mode energy observed during the sequence, 
which is assumed to be similar to the energy at the source. Only mode 1 is presented 
here. Flowever, the results for mode 2 and 3 looked quite similar. 

Figures 14, 16 and 18 show the 400Hz component of the mode energy for the 
three case studies. Instead of working with the whole broadband pulse, this study chose to 
extract the 400Hz component from the signal for the initial analysis; the narrowband 
comparison is then simpler since the spreading induced by the non linear coupling 
between modes is removed. Additionally, this central frequency carries the main acoustic 
information. 
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Figure 13. Relative Mode energy for mode l(in dB) received at the VLA between 
1815 GMT on May 5, 2001 and 0145 GMT on May 6. A few soliton locations are 

indicated with the red lines. 
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Finally, this section compares the results obtained during the first and second case 
studies, when solitons were present in the acoustic path, to the results observed during the 
third case, which this work has called quiet time. 


^mpfituec squared of node 1.2 and 3 measured atif>c VL*\ ^ frequertoy oomiMoent: ^OOHz 





Figure 14. Temporal variations of the 400Flz component of the mode acoustic energy 
recorded at the VLA for mode 1, 2 and 3 between May 5, 2001 at 1945GMT and 

May 6, 2001 at0145GMT. 
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Starting with this first case study, a few main characteristies are readily 
notieeable. First, as shown in Figure 13 eompared to Figure 17, the signal seems to be 
more spread in time. Moreover, the analysis of the 400Hz signal (Figure 14) emphasizes 
the faet that the signal fiuetuates more quiekly as an ISW enters the aeoustie path. Then, 
in the presenee of solitons, mode energy also seems to deerease on average, as illustrated 
by the eomparison of Figures 14 and 18. On May 5 (Figure 14) and 6 (Figure 16) mode 1 
and 2 seem to convey a little bit less energy in the presenee of solitons than during the 
quiet time (Figure 18), while mode 3 seems to earry more energy. Many of the referenced 
papers talk about mode eoupling and deseribe how this phenomenon was observed while 
dealing with aeousties and solitary internal waves; this theory also forms the eore of the 
model presented in Chapter II and is responsible for the exehanges of energy between 
modes. Rousseff et al. (2002) detailed these exehanges of energy, but this paper foeuses 
only on the predictability of these exehanges rather than quantifying the observed 


eoupling. 
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Figure 15. Same as Figure 13 for a signal observed between 0700 GMT and 

1500GMT on May 6, 2001. 


32 





















Ampfetude squared of mode 1.2 and 3 measured art the VLA — frequence compor>ent: 40OI-tz 



Figure 16. Same as Figure 14 for a signal observed between 0700 GMT and 

1500GMT on May 6, 2001. 

Besides the exchanges of energy between modes, the variability of the signal also 
signifies the presence of solitons along the acoustic path. In fact, some evidence of 
variability is clearly observable on May 7, during the so-called quiet time, though the 
fluctuations measured on May 5 and May 6 are much faster. Indeed, as stated in chapter 
1, non linear internal waves are not the only parameters which control the variability of 
the water column. In contrast, linear internal waves, as well as subtidal or tidal 
phenomena, occur simultaneously and are probably responsible for the changes in the 
mode energy for the third case study, on May 7. 
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The comparison between May 5 and May 6 also reveals some interesting features. 
Surprisingly, the presence of more than one soliton on the acoustical path does not seem 
to significantly affect the amplitude of the mode energy oscillations observed at the VLA. 
Compared with the soliton on May 5, the three ISWs on May 6 were also bigger. Based 
on the theoretical model presented in the predictions outlined in Chapter II, one could 
expect more dramatic effects in the presence of these solitons. However, it is probable 
that the presence of three solitons makes the constructive and destructive phases more 
complicated and may explain the previous result. 
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Figure 17. Same as Figure 13 for a signal observed between 0315 GMT and 

0715GMT on May 7, 2001. 


Finally, from these few pictures it appears that non linear internal waves might be 
partly responsible for the time spreading of the signal. This theoretical model predicted a 
small amount of spreading caused by solitons, and most of the spreading could be 
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attributed to the presence of other phenomena, such as linear waves or mesoscale events; 
this point will be further developed later in this chapter. 





Figure 18. Same as Figure 14 for a signal observed between 0315 GMT and 

0715GMT on May 7, 2001. 

B. COMPARISON BETWEEN MODEL AND OBSERVATIONS: MODAL 
ENERGY VARIABILITY 

As stated in the introductory chapter of this report, the goal of this project was to 
test the fast theoretical model designed by Colosi against the in-situ data presented in the 
previous section. The following results were evaluated from MATLAB based programs 
built upon this theory. The runs were conducted both on UNIX and Windows systems. 
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Before looking elosely at the results, this seetion eonsiders a few diffieulties 
regarding the strengths and weaknesses of the eomparison. First, as explained in ehapter 
II, this model was designed for quiek implementation, and many assumptions were made 
for this purpose, obviously deereasing the aeeuraey of the study’s results when the mode 
eoupling is strong. Seeondly, as diseussed earlier, this model negleets tidal and subtidal 
variabilities, as well as the effeets eaused by random linear internal waves, though these 
elements eould signifieantly alter the observed aeoustie signal (Chiu et al., 2004). 
Thirdly, data eolleetion methods applied during this eampaign of measurements did not 
perfeetly mateh those required to exhaustively evaluate the model. For instanee, the 
loeations of the different oeeanographie moorings made a few parameters, partieularly 
ISW speed, very diffieult to evaluate. This parameter is, however, essential in the 
evaluation of the eharaeteristie half-width of the soliton and is a sensitive input for this 
model, as will be demonstrated later. Finally, bathymetrie eoupling ean also influenee the 
ineidental mode amplitudes on the NLIWs and thus represent an additional eomplieation. 

For any of the three test ease studies, this study eondueted a direet eomparison 
between the data and the model as illustrated in Figures 19, 21, 23, 25 and 27. In these 
figures, the 400Hz eomponent of the mode energy reeeived at the VLA is displayed with 
respeet to the loeation of the soliton in the aeoustie path. The modeled data in these 
figures are represented with blue lines and blue dots, while the observed data are in red. 
Signifieantly, in these five figures, like the observations, the model demonstrated fast 
ehanges in the predieted energy reeeived at the array. Although these results were 
promising, a direet superposition of the observed data on top of the modeled data 
appeared diffieult. Indeed, many of the statistieal tools employed here are hypersensitive 
to slight ehanges in the inputs, sueh as the soliton speed, width or amplitude. Therefore, 
this study superposed only a eouple of signifieant points over the simulated data for any 
sequenee of 88 pulses. In addition to these points, whieh represent averages over the first 
and seeond parts of eaeh sequenee, errorbars eneompass the extremes on both sides of the 
mean values. Finally, all the observed data were sealed with estimates of the unperturbed 
mode amplitude, evaluated from seetions where there were no NLIWs. Figures 20, 22, 
24, 26 and 28 are seatterplots that eombine observed and modeled data. The value of the 
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correlation coefficient is indicated in the upper part of these figures in order to help us in 
assessing the quality of this model and the relevance of the comparison. These latter 
values were also compared to the sample error s defined by ^« -silln, where n 
represents the number of events accounted in the correlation calculation. Finally, the 
soliton amplitude, speed and width used as inputs for the runs are similar to the 
information presented in Chapter 11. 


May 9tH 




Figure 19. Relative mode energy fluctuation conveyed by the 400Hz carrier 
frequency and received at the array for mode 1 and mode 2 on May 5, 2001. The 
blue lines, dots and error bars represent the modeled data while the red dots and 

whiskers are the observed ones. 
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Figure 20. Scatterplot comparing observed mode energy to modeled mode energy for 
mode 1 and 2 on May 5, 2001 for the 400Hz carrier frequency. 

A few main results can be extracted from the figures above. First, as expected 
from the definition of the mode amplitude described by Equation 2.17, the simulation 
predicts huge oscillations of the mode energy. Flowever, because of the changes in the 
bathymetry and in the shape of the soliton, these oscillations lack periodic consistency 
and show some variability as the soliton moves shoreward. While the model successfully 
predicts the mode energy variability for mode 1, its results are quite poor regarding mode 
2 . 

Obviously, the results are not perfect, and the correlation coefficient is just a little 
better than the sample error, ^ »0.35. Nevertheless, as illustrated in the scatterplot, 
almost all the points are in the neighborhood of the slope 1 line, and the correlation 
coefficient is decreased by a few outliers. For mode 2, it seems difficult for the model to 
predict the right energy, as illustrated in Figure 20. Figure 19 suggests that the model is 
able to track the tendency of the energy changes, but the continual presence of time leads 
and lags makes the statistical results quite poor. Many factors could explain these 
differences between mode 1 and higher order modes, particularly the fact that this model 
uses a weak, single scattering approximation which may no longer hold for modes higher 
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than one. Moreover, random linear internal waves, the eomplex bathymetry or other 
kinds of variability can induce mode coupling, and are thought to affect these higher 
modes more. 

Finally, the model appears to generate better results in the first part of the plot 
when the soliton has not yet reached the sharp continental shelf between S5 and S4. This 
result could have been anticipated, since the model does not take into account the non 
linear interactions of the soliton with the bottom. 


May 5«h 




Figure 21. Same as Figure 19 for the 395Hz carrier frequency. 


This study gauged the ability of this model to predict not only the energy 
conveyed by the main carrier frequency, but any component of the lOOHz-bandwidth 
pulse. To this end, this study conducted the same kind of analysis for two other 
frequencies, namely 395Hz and 405Hz, which also carry a lot of energy. Figures 21, 22, 
23 and 24 present these results, which bear similarities to those observed for 400Flz. The 
fourth section of this chapter explores the sensitivity of the model regarding frequency. 
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Figure 22. Same as Figure 20 for the 395Hz earrier frequency. 
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Figure 23. Same as Figure 19 for the 405Hz carrier frequency. 
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Figure 24. Same as Figure 20 for the 405FIz carrier frequency. 

After testing the model’s energy predictions, the same kind of analysis was 
conducted for May 6 as a series of three solitons propagated along the acoustic path. 
Technical failures prevented proper collection of acoustic data during the whole period of 
time, so only a few points were available for this analysis. Flowever, in this case, which 
used only ten points to evaluate the results, the correlation coefficients are even higher. 
Then, the expected value for the sample error is given by ^ « yj2/n = 0.45 . These values 
confidently suggest that the model’s results are relevant, and that there is a significant 
difference between results between modes 1 and 2. Moreover, the fact that the coefficient 
correlation reached a value close to 0.8 in the second case study is very promising for the 
model. Because three large solitons were observed along the acoustic path in this case, it 
is likely that their influence on the water column variability is much bigger than for the 
first case study and largely overcomes the relative influence of other environmental 
variability. 
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Figure 25. Same as Figure 19 on May 6, 2001. A series of three huge solitons was 

observed along the aeoustie path. 




Figure 26. Same as Figure 20 on May 6, 2001. 
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c. 


COMPARISON BETWEEN MODEL AND OBSERVATIONS: TIME 
SPREADING 


Part of the effects caused by the solitons on the signals observed at the array was 
related to energy fluctuations arising from mode coupling, as demonstrated earlier. 
Looking at Figures 13, 15 and 17, it is nevertheless notable that the signal shows 
significant time spread in the presence of solitons relative to the case without ISWs. In 
fact, mode coupling is known to act on modal phase as well as modal amplitudes. The 
current section focuses on this aspect of the interactions between solitary internal waves 
and sound propagation. Then, the time spreading observed in the in-situ data were 
compared against the simulated results. Two main variables were used to estimate the 
spreading; the time lagged-intensity covariance of the signal and the travel time 
fluctuation. The first variable more adequately characterizes the spreading, which affects 
in-situ data, and is described by the following equation; 


= [ + where 4(r) = |A„(r)| 


Regarding the spreading predicted by the model, chapter II demonstrates that the 
mode amplitude phase fluctuation is frequency dependent and given by; 


k 

j=\ m=l 




Then, the travel time fluctuation of the mode is given by; 

^_diS0„{R)) 


If the phase varies linearly across the bandwidth of the signal, the pulse does 
not spread, but is only advanced or delayed. This advance or delay is termed a travel time 
bias, Tj. 
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Figure 27. Simulated phase aeross the 350 to 450 Hz bandwidth on May 5, 2001. The 
soliton is located 30km away from the source. The linear fit to the phase is shown 

in green. 


However, if the phase variation is non-linear, then dispersion or time spreading 
occurs. Figure 27 shows an example of the phase variation across the frequency band 350 
to 450Hz for the simulations on May 5; clearly, both modal bias and spreading are 
occurring. The modal time bias is qualified by fitting a line through the phases <76*,^, 


and the bias thus becomes the slope of this linear fit. To quantify the spread, one subtracts 
from the linear fit, forming a detrended phase function 50^, which is then used in 

Equation 3.4 to compute time fluctuation as a function of frequency. The pulse spread, 
which is also the root mean square variation of r across the bandwidth, can be expressed 
by: 


T 


0 



(3.5) 


Figure 28 shows the estimates of the bias and spread as a function of the ISW 
position Tq . What is first evident is that the bias and spread increase when the non linear 
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wave is more distant from the souree. This is understandable, since the modes propagate 
at different speeds, and thus accumulate larger time differences at larger ranges. The time 
fluctuation of a mode will be greatest if the observed energy in that mode travels as much 
of the source receiver distance as other modes and only couples near the receiver. The 
modulation of the bias by the beat wavenumber resonance is also notable. Depending on 
the relative phase of the modes at the ISW location r^, the bias can be either positive or 
negative. Additionally, the bias and spread represent significant fractions of the 400Hz 
period (2.5 ms), so that these corrections could be very important for time coherent 
processing. Further remarkable are results from simulated data based on soliton 
characteristics from May 6, as shown in Figures 29 and 30. Clearly, these last two 
figures, particularly the values of the spread, suggest that the signal would experience 
less spread as a result of many large solitons rather than a single isolated one. 



Figure 28. Mode 1 pulse bias and spread on May 5 as a function of the distance from 

the source, . 
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Figure 29. Same as Figure 27 for May 6. 



Figure 30. Same as Figure 27 for May 6. 
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As illustrated in Figure 31, and as predicted by the model, the observed spreading 
is quite noticeable for the case studies and more dramatic on May 5 than on May 6. The 
model also predicted a larger spread of the signal as the NLIW neared the array. This 
result was unfortunately not truly observable from the data, as illustrated in Figures 32 
and 33. Finally, it is significant that this spreading caused by the ISW is relatively small, 
less than 3 ms, and does not explain all the distortion in Figures 13 and 15. 

Figure 31 displays the time lagged-intensity covariance of the signal observed for 
the three case studies and focuses on the comparison between these three days for mode 1 
and mode 2. Figures 32 and 33 more carefully detail the evolution of the spreading with 
time. Indeed, in these figures, the case study is separated in two halves (1®* part and 2"‘* 
part) to track the evolutions caused by the propagation of the soliton. The first part (red 
curves) represents the energy observed when the soliton was more than 15 km from the 
array; the second part (green curves) characterizes the period when the soliton was less 
than 15 km from the array. Finally, the variables cTq, cTjj and CTjg characterize the e- 
folding of the covariance for the three case studies and are known to characterize the 
spreading. In Table 1, the spreadings observed on May 5 and 6 are represented by 
%5 and , defined by: 


a\ 


sp5 


*^15 ■‘^0 


and CTjpg CTjg-CTQ 


Table 1. Covariance cTq, and cr^^g. 



May 7, 2001 - cTq 

May 5, 2001 - 

May 6, 2001 - cr^^g 

Mode 1 

6.1 ms 

10.7 ms 

6.1 ms 

Mode 2 

6.4 ms 

10.5 ms 

5.0 ms 

Mode 3 

6.5 ms 

10.2 ms 

3.9 ms 


47 






Figure 31. Time lagged-intensity covariance of the signal on May 5, 6 and 7, 2001. 

As evident in Figure 33, the analysis was unable to observe any difference 
between the first part of the path and the second part, as expected from the modeled 
spreading, despite the presence of a soliton in the acoustic path. Moreover, the observed 
spreading was definitely larger than the modeled. From these observations, spreading 
caused by the ISW appears to be much less important than the spreading caused by other 
environmental variability, particularly the random linear internal waves. 
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Figure 32. Time lagged-intensity eovarianee of the signal on May 7, 2001 for mode 1 

and mode 2. 




Figure 33. Same as Figure 32 on May 5, 2001. 
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D. SENSITIVITY OF THE MODEL 


The last step in the evaluation of the model was to test its sensitivity against some 
of its inputs, namely the propagation speed of the soliton, its amplitude and width and a 
few frequencies within the bandwidth. 

As the environmental conditions, such as the temperature, salinity and density, 
were quite similar during the setting of the solitons on May 5 and 6, 2001, the model was 
able to investigate the relative influence of the number of solitons, their amplitude and 
width by just looking at the simulation for these two days. The propagation speed was 
also assumed to be identical in the following comparison. 





Figure 34. Simulated modal energy received at the VLA on May 5 and 6, 2001. 
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Clearly, any of the three solitons observed on May 6, 2001, was bigger than the 
soliton reported on May 5, 2001. In fact, the amplitude of the soliton does not change the 
whole shape of the energy variation; in contrast, the amplitude tends to magnify the 
existing variations, as illustrated in Figure 34. The addition of multiple waves tends to 
create some constructive and destructive effects in the simulation, which may explain the 
observed lags and leads between the two simulations. 



Figure 35. Five different ISW speed profiles between moorings S7 and S5. 


This section also examined the propagation velocity of the soliton. Actually, 
during the numerous tests conducted for this project, the model/data comparison was 
highly sensitive to this variable. The ISW speed was also a critical factor since it was 
difficult to pinpoint accurately. Indeed, the only available and appropriate data were the 
currents at moorings S4, S5 and S7, the temperature measurements at these same 
locations, which revealed the signature of the soliton, and a few SAR images, which 
indicated the direction of propagation itself From these different observations, the actual 
speed of the soliton could be deduced, followed by the speed of the soliton as it is seen 
across the acoustic path. This latter speed proved to be slightly bigger since the 
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propagation of the soliton and the acoustic path had different orientations. Then, different 
speeds were simulated, as illustrated in Figure 35, and the model was run for any of these 
five simulations. 


May 5«h 




Figure 36. Simulated fluctuations of the 400Hz component of the mode energy for 
mode 1 (upper panel) and mode 2 (lower panel) as a function of the propagation 
of an ISW in the acoustic path for 5 different ISW speed profiles. 


These examples were used to test three linearly decreasing speed profiles 
associated with the motion of the soliton from S7 to S5, two possessing gentle slopes and 
one with a steep slope. This work also simulated a constant speed profile and a depth- 
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dependent profile that agreed with the linear wave theory in shallow waters. As 
demonstrated in Figure 36, the differenees between the five runs are dramatie, and both 
the average speed over the path and the slope seem to be important. 
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Figure 37. Simulated fiuetuations of the mode energy for mode 1 (upper panel) and 
mode 2 (lower panel) as a funetion of the propagation of an ISW in the aeoustic 
path for three different frequency components of the signal; 375Hz, 400Hz and 

425Hz. 

The last variable in the critical analysis of the model was its frequency 
dependence. As illustrated in Figure 37, and earlier in Figures 21 to 24, the results were 
conclusive. Furthermore, the model itself proved greatly promising, as it was able to 
predict the energy variations for any frequency component of the signal. 

E. DISCUSSION AND COMMENTS 

Despite the difficulties in extracting the appropriate data for the 2001 Asian 
experiment, the model proved highly adept at predicting the energy variations when a few 
solitons came into the acoustic path. This was particularly true for mode 1, which was the 

easiest one to resolve with the available vertical array. The results observed for mode 2 
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and for any higher mode did not meet the study’s expectations; as stated earlier, many 
explanations could be offered to explain this unforeseen result. For one, coupling effects 
and scattering effects are much more complex for these higher modes, which are also 
probably more sensitive to fluctuations caused by linear internal waves, the bathymetry 
or other environmental variability. 

Regarding the time spreading, the modeled data did not agree with the observed 
data, and once again these differences are most likely the result of influences by other 
kinds of environmental variability. 

Finally, despite the vast amount of data collected during ASIAEX 2001, this study 
did not collect adequate ISW speed profdes; unfortunately, this data is essential to run the 
model. Thus, with more constraining data, the results from the model could have been 
different. Future runs of the same kind of comparison should consider working a shorter 
range for the acoustic path and a simpler bathymetry, which would suppress some 
undesirable variability and coupling. 
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IV. DETECTION OE NON LINEAR INTERNAL WAVES BY 
SYNTHETIC APERTURE RADARS 


Internal waves are one of the first oeean phenomena detected on synthetic 
aperture radar (SAR) images, where they appear as a sequence of bright and dark quasi- 
parallel lines. This result was early described by Alpers (1985), Apel et al. (1979). Based 
on the quality of these pictures, very promising results are consistently observed from 
SAR satellites. In order to understand the strength of these systems in characterizing 
solitons, this next section provides a review of some central concepts concerning SAR 
images. Following that is a review of the basic theory supporting the detection of solitons 
by SAR systems. Finally, the study presents an elegant theory by Zheng et al. (2001) that 
proposed two simple ways of retrieving the characteristic width and amplitude of the non 
linear internal waves. The following chapter determines the relevance of the theory by 
testing it against in-situ data for different environments. This analysis gives particular 
attention to wind effects. 

A. SYNTHETIC APERTURE RADARS AND OCEANIC PURPOSES 

Synthetic aperture radars were developed in the 1970s in order to improve the 
spatial resolution offers by real aperture radars. The first satellite to carry a SAR system 
was SEAS AT (1978), which provided very promising results regarding wind retrieval, 
detection of oceanic features such as oil slicks and non linear internal waves and ship 
tracking. These abilities are particularly important because clouds are transparent for 
radar satellites and radar information can be retrieved either by day or by night. 
Moreover, radars are active systems and allow advanced data processing for finer 
resolution. Finally, as their temporal and spatial coverage and storage capacities improve 
in the future (RADARSAT 2 vs. RADARSAT 1), they will become more and more 
applicable. 

From real aperture radar physics, it is known that for a specific frequency / (or 
wavelength X ) and a slant range R , the azimuth resolution is entirely dependent on the 
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aperture length . However, for the altitude at whieh satellite imaging sensors operate, 
engineering difficulties hamper achievement of L^l ^ values greater than several 

hundreds, this value being inversely proportional to the azimuth resolution. Through the 
use of synthetic aperture radars, however, 10 meter-resolutions can be achieved. In the 
1950s, it was realized that the Doppler spread of the echo signal could be used to 
synthesize a much longer aperture so as to greatly improve the resolution of side-looking 
radars. It is beyond the scope of this chapter to consider all the applications and 
theoretical concerns behind SAR images; this chapter will focus only on the theory that 
supports soliton detection by SAR images. Table 2 summarizes the characteristics of the 
main SAR systems. 

Table 2. Characteristics of the main SAR systems. After Olmsted (1993). 



y.fKsion 

Orbit 

EKS-1 

JERS-1 

Radaisat 

Alfinir^ 


sft,S km 

/9 3-.S ’I krn 

lucluutiou 

98.516“ 

97662“ 

98.594“ 

Pciiod 

<5327 907 s 

5799-72 s 

<5042 s 

Rjcpear C>xle 

3. (35. 176) davs 

44 davs 

24 davs 

Orbit 

sun svnclironous 

sun s\Tichronous 

sun ssTichronous 

Ground Track V^elocit>' 

6 628 km/s 

6 883 km s 

6.576 km s 

Instrument 

l-reqti«ic>' 

t. -Band O i GHz) 

L-Band (1.2/5 GHz) 

C-Btnd (5.3 GHz) 

Wavelen^ 

5.66 cm 

23 5 cm 

5 66 cm 

Pulse Repeitiou Freq. 

1640-1720 Hz 

1505 8-1606 Hz 

1270-1390 Hz 

Pulse Leuflth (BW;) 

37.1 MS 

35 MS 

42 MS 


(15.5 MHz) 

(15 MHz) 

(11.617.3J0 MHz) 

Polarization 

W 

HH 

HH 

.Antenna SizeL x W 

lOX Im 

11 9X2 4m 

15 X 1 5 m 

PeaK Power 

4 .S k\V 

1.5 kW 

5kW 

Average Power 

300 \V 

71 W 

300 R* 

Noise Equivalent O’ 

-18dB 

-20.5 dB 

-21 dB 


Swatli VVkUJi 




Max Resolution Rg x Az 



lOX 10m 

Resolution (3 # looks 



23X30t? 4 

^^•sTern 

On Board Storage 

None 

Tape 20 mm 

2 Tapes ^ L4.5 mm 

Look Angle 

RiiOil 20.355* 

RmUi 35.21* 

R & L 20-50* 

Incidence Allies. Mid 

19.35-26.50“. 23“ 

36.14^1 51“, 38.91“ 

22.64-59.56“, 45.12“ 

kootpnnt K.f X Az 

8t'X4.S5km 

/OX 14 km 

50-150 X 4.3 km 

Doppler BW 

1260 Hz 

1157 Hz 

939 Hz 

Coherent Integ Time 

.6 % 

1 7 s 

46 s 

Windc-wmc 

Hamming 

Hamming 

Hamminc 

Pulse Compression Ratio 

580; 1 

525;1 

491. 727. 1260; 1 

Range Sample Rate 

13.96 MHz 

17.1 MHz 

12.9, IS.5, 32.3 MHz 

Maxuuiuu Data Rate 

105 Mb s 

30 Mb X 2 Mb s 

105 Mb s 

Quantizatioii 

5 bits sample 

3 bits sample 

4 bits sample 
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For a homogenous medium with a slightly rough surfaee and an intermediate 
ineident angle, generally between 20° and 70°, the seattering from the surfaee ean be 
deseribed using the Bragg model. The Bragg seattering theory and its domain of 
applieability are deseribed eomprehensively in Plant (1990), and only the most important 
results for this projeet will be presented here. This theory states that any seetion of that 
surfaee wave speetrum resonates with the ineident wave in produeing strong 

2L 

baekseattering. The Bragg resonanee eondition is also given by: —.sin 0 = n where L is 

X 

the oeean wavelength, X is the radar wavelength and 0 represents the ineident angle 
assoeiated with the satellite system and n = 1 , 2 ,... 

Radar images are also eonsidered to be two-dimensional representations of the 
intensity of radar return signals. If one assumes that the radar reeeiving system is a 
perfeetly linear system, the intensity of radar return signals should depend linearly on the 
roughness of the target surfaee, whieh is quantified by the baekseatter eross seetion per 
unit area. For an oeean surfaee. Plant derived an expression for the normalized radar 
eross seetion of a slightly rough surfaee, cTq : 

=16.;r.k^|g,(^)f .^2^08111^0), (4-1) 

where kg is the wavenumber of the eleetromagnetie signal. ^^( 0 , 2^0 sin6()) is the two- 
dimensional varianee speetrum of surfaee displaeement evaluated for the waveveetor 
{k^ = 0,ky = 2ko.sin65,). The axes of referenee are relative to the viewing direetion of the 

satellite, the x-axis being parallel to the traek of the satellite. In other words, the 
baekseatter is eaused only by surfaee waves whieh lie within a restrieted wavelength 
range and travel nearly toward or away from the antenna. It is further notable here that 
this expression breaks for steep waves. For this latter ease. Plant (1990) proposed a mueh 
more eomplieated equation ealled the eomposite surfaee theory, whieh ineludes the mean 
surfaee slope and displaeement. Finally, the indiees ij denote the polarizations of the 

ineident and baekseattered radiation, and, if is the relative dieleetrie eonstant of 
seawater, gyiO) is the first order seattering eoeffieients given by: 
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for horizontal polarization and, 


(4.2) 



Thus, this set of equations states that, for a given radar wave number 
incident angle 6 ^^, the intensity of the radar return signals depends only on the two- 

dimensional wavenumber spectral density of the surface ocean wave field, 'F. As 
indicated in Table 1, the Canadian SAR system RADARSAT, which is used in the 
analysis, is polarized horizontally. 


B. INTERACTION BETWEEN NON LINEAR INTERNAL WAVES AND 
SURFACE WIND WAVES. 


Based on measurements by Zhang and Cox (1994), and the balance equation for 
the two-dimensional wavenumber spectral density of the surface ocean wave field, Zheng 
et al. (1994) proved that 'F could be written; 



Even if this equation is quite coarse, the method is still elegant. Indeed, Zheng et 
al. (2001) retrieved 'F from the principal of conservation of wave energy. Actually, the 
balance equation for the two-dimensional variance density spectrum of surface elevation 
can be written; 

^ + (y + c, ).V„>P = S., + S,„ + S„ + . (4.4) 

In this equation, largely used in wave models, 'F is calculated for any frequency 
/ or wavenumber k and any direction (p. and represent the source 

terms for wind wave interaction, resonant wave-wave interactions, dissipation and wave 
current interactions. Assuming a steady state and using previously derived solutions for 
the source terms, Zheng et ah, retrieved equation 4.3. 
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The coefficient originates from the wave-wave non linear term presented by 

McGoldrick (1965); this project does not directly use the coefficient and therefore will 
undertake no further investigation of it. The value of m was equal to 0.04, m * is the 
friction velocity, y is the viscosity of seawater co is the dominant angular frequency of 
the surface waves. S^p is the excess momentum flux tensor, [/^is the velocity 

component of the large scale current field , k is the observed wavenumber. The third 
term on the right side of equation (4.3) is given by; 


’ap- 


dx„ 


1 

2 


.du 2 , dv. Sv . 2 

(—.cos f + {—H-).cos^z).sm^z)H-.sm (p) 

dx dy dx dy 


(4.5) 


u and V were said to be the velocity components of the surface currents and (j) is the 
wave direction, a and (5 are subscripts representing the horizontal axis of reference. 
This latter term, illustrated in Figure 38 and Figure 39, relates to the objectives of this 
theory and will be investigated more closely in the next section. 

The manipulation of the previous equations requires careful handling of the 
different coordinate systems used by Zheng et al. Plant derived equation 4.1 using a 
coordinate system related to the viewing direction; equation 4.4 must be used in a 
Cartesian coordinate system associated with the cardinal directions. Finally, Zheng et ah, 
used equation 4.5 in a coordinate system associated with the direction of the non-linear 
internal waves. Despite the confusion in the notations, the main ideas behind this theory 
offer valuable qualitative estimates of the influence of the different components on the 
backscattered energy. The section dedicated to the wind influences looks more closely at 
the definition of . 


Because the spatial variability of the viscosity is small in the oceans (except in 
estuaries, coastal regions, polluted areas, and frontal zones), the contribution of the 
second term in equation 4.3 to a radar image can be considered as a constant. Obviously, 
the two other terms play key roles in generating an ocean radar image. 
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Figure 38. Schematic of the interaction between non linear waves and surface wind 
waves and their effects on SAR images. After Hsu and Liu (2000). 


It is significant here that equation 4.3 is based on the principal of conservation of 
wave energy and does not hold for rapidly changing wind conditions. Nor does it hold in 
very shallow waters, where dissipation or non linear interactions can be very important. 
Moreover, this equation was designed for gravity-capillary waves and a different 
equation should be used for capillary waves. However, regarding the wavelength of the 
most common SAR systems, such as ERS and RADARSAT, this equation is sufficiently 
accurate. 


Surface soliton 


R i p 


Internal soliton 



\V\\\ \ \\ \\\\\VA \\'\ \\\\ \ \ \ \ 


Figure 39. Internal soliton in a two-layer fluid with finite depth. Dashed lines are 
isotachs of the water particle speed. Arrows indicate magnitude and direction of 
the water particle movement. After Osborne and Burch (1980). 
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C. DETERMINATION OF THE CHARACTERISTIC HALF WIDTH AND 
THE ASSOCIATED AMPLITUDE IN THE KORTEWEG AND DEVRIES 
MODEL 


In their paper, Zheng et al. (2005) presented two methods based on the KdV 
equation. This section also summarizes their basic ideas, while the following chapter 
reports the results that were found using their techniques. 

The first chapter of this report introduced the KdV equation for a two-layer 
model, and its solution was described by a secant function, such that; 

?l{x,t) = %sech\^^^) (4.6) 

A 


where A = 



is the characteristic width and V =Cq + 


«-?7o 

3 


the soliton velocity. 


The first technique presented by Zheng et al., is based on the direct calculation of 
the interaction between the soliton and the surface waves. Then, using the previous 
definition of the depression in a barotropic two-layer model, one can deduce the velocity 
in the upper and lower layers, respectively Wj and Wj, such that; 


Wi(v,0 = (-V^).sech ( —-—), 

/ij A 

Wj (x, t) = -( ° ° ). sec h (-). 


h 


A 


(4.7) 


One can then write that; 


^ B Cn.?7n.COS^ (zi 

= —— - -. tanh( 


'ap- 


dx„ 


1\A 


A A 


(4.8) 


Using equation 4.1 and 4.8, one derives an expression for the backscattered 
energy received by a satellite and caused by a soliton; 


{0),j=Q{k,gi,,9,(pA)-[ 


C,. 71 ,.COS^ (!) , , 


,x-V.t, 


. tanh(^;^-^). sec ] 

fi^.A A A 


(4.9) 
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with Q{k,g.j,O,(p,d) = \6.n.kl.\g^j{0i\ .f{(p,d). 


(4.10) 


Essentially, this equation means that it is the alternating compressive and tensile 
effects of the internal wave currents that render them visible on the surface. Actually, the 
modification of surface wave spectrum by internal waves can be very significant, 
increasing the spectral level on the order of 200% (Apel, 2007). Moreover, for a SAR 
satellite, the frequency and the antenna polarization are fixed, and the incidence angle, 
the viewing direction, the propagation direction of the soliton and the dielectric constant 
change only a little within the range of an ocean internal wave packet. Therefore, Q can 
be treated approximately as a constant over the range of a soliton. Equation 4.10 is 
slightly different from that proposed by Zheng et ah, and takes into account the looking 
direction of the satellite as well as the propagation direction of the non linear internal 
waves through the term f{(p,d). If the latter need not be thoroughly evaluated in this 
theory, it is, however, essential to cover while looking at the wind effects. 

In order to retrieve the characteristic half-width, one can use a normalized version 
of equation 4.9 and compare it to the normalized and centered SAR pixel values along a 
transect profile. A curve fitting method would also render the raw data smoother. 
Because SAR images inherently contain noise, this work relied on hand chosen pixels 
that characterize the soliton. They roughly represent the major features of the bright and 
dark lines associated with the soliton. 


As the next chapter discusses, this method offers interesting results that can be 
altered by wind speed and direction. To overcome this difficulty, Zheng et ah, introduced 
another method, which they called the peak-to-peak method. This method aims at 
retrieving the characteristic half-width by measuring the distance between successive 
maximum and minimum. Equation 4.9 proves the existence of these remarkable points, 
which are characterized by bright and dark lines on SAR images. Mathematically, these 
two points satisfy the following equation; 


dcT, 


o./vr 


dx 


= 0 


:tmax 

xmin 


(4.11) 
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which is equivalent to; 


ry pp • ry pp 

SQch\^ - )-ltgh\^ -) = 0 (4.12) 

A A 

with being the distance between the maximum and the minimum. 

Solving this equation, one obtains the following relationship between and A; 

</„=2ln(i^^).A = 1.32A. (4.13) 

This demonstration, derived by Zheng et ah, forms an essential foundation for 
understanding the whole theory of this work. Because it requests only two points, this 
method is very powerful and overcome most of the diffieulties related to other factors 
such as wind. As later sections will discuss, however, the minimum is diffieult to identify 
in the presence of strong winds or small ineidenee angles. 

From these two methods and the KdV theory presented in the first chapter it is 
then possible to retrieve the amplitude of the perturbation based on the definition of the 
characteristic half-width given by equation 4.6. 

D. COMMENTS ON WIND EFFECTS 

Equation 4.3 provides an elegant means to qualitatively estimate the influenee of 
different variables on the radar cross-seetion measured by SAR systems. As explained 
above, this work simplifies many of the previous equations and omits several important 
features of the wind effects. In fact, Zheng et ah, limited their analysis of the wind effects 
to comparisons of orders of magnitude between the wave-eurrent source term and wind 
input source terms. 

Wright et al. (1977) first introduced the expression for the wind input presented 
by Zheng et al. Sinee then, Chaltkov and Belevich (1993) have designed a eonsiderably 
more complicated expression, currently in use by advanced wave models sueh as WAM; 

= (4.14) 
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F(/,^)is the two-dimensional varianee density speetrum of surfaee elevation 
expressed as a funetion of wave frequeney and direetion. This funetion is equivalent to 
the previous 'F . (5 is a non dimensional wind-wave interaetion parameter whieh ean be 
approximated by; 


Q, =1.075 + 75Q 
Q2=1.2 + 300Q 

Uq = 0 . 25(35 ! 


\Q\/3 = 




ayO)^.{a,.o)^-a,)-a^ 

for 

(a,.S)^-a^) 


<f2j/2 

Qj / 2 < < Qj where 

Qj < < Qj 


ai=0.25 + 395Q 
a2=0.35 + 150Q 

Uj =((3o -(32 -(3j)/((3o -(34 +( 35 ) 


a4=0.3 + 300C2 

(35 =(34iTj 
Ug =( 3 Q.(l-a 3 ) 

= (flg-CQj “1)^ + ‘ 2 io )/(^2 

(3g =(32.Q, 


(4.15) 


^9=0.35 + 240^ 
Oio =-0.05 + 470C2 


In this set of equations, o)^ = cos{(p-(p^), it is the non dimensional frequeney of a 

speetral eomponent . (p^ is the wind direetion, is the wind veloeity at some partieular 
In 


elevation X = 


fc|oos(^-^„) 


and finally, is the drag eoeffieient. 


This system is obviously quite eomplex but presents the advantage of assessing 
the wind influenee on baekseattered signal through the term (5. As illustrated earlier, the 
return radar signal depend on the viewing direetion from the satellite (p and its 

^n 

wavenumber k = 2.kQ.sm0 = —.sin6*, with 6 being the ineidenee angle of the satellite 

Xq 

system. As the frequeney is fixed for a partieular radar system, the wind ean affeet the 
results beeause of its magnitude and the differenee between the looking angle cp and 

the wind direetion Thus strong wind and upwind/downwind looking angles will favor 
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a strong background radar signal. This effect can obviously cover the signature of non 
linear internal waves or prevent diserimination of minima. In their paper, Zheng et ah, 
assessed that the maximum wind at whieh soliton deteetion by SAR systems is possible 
should be around 5 m/s. Aetually, the maximum wind is thought to be higher than that, 
around 25m/s for ERS-1 (Raney, 1995) and around 15 m/s for RADARSAT (Staples, 
personal eommunieation, 2007). This limit depends on the soliton direetion and strength, 
the wind speed and direetion; it is also highly dependent on the speeifieations of the 
satellite system, partieularly its level of saturation. 

A lower threshold eould be defined as well. Indeed, as presented in the Bragg 
model, the deteetion of oeeanie features by SAR systems is based on their analysis of the 
wave speetrum for the radar wavelength. Thus, a smooth oeean observed during no wind 
eonditions will not baekseatter anything toward the satellite system, but will refleet the 
radar signal away from it at an angle equal to the ineidenee angle, as stated by the laws of 
opties. The wind must also be strong enough to develop eapillary waves or small gravity 
waves deteetable by SAR satellite. From wave theory and observations, this limit eould 
be reaehed for wind speeds greater than 1.5 m/s. Additionally, a seeond parameter 
eontributes to the definition of this lowest limit, namely the deteetion threshold of the 
satellite system. This parameter, whieh is highly dependent on the satellite, proves 
diffieult to evaluate. In faet, as stated by Raney (1995), the minimum wind speed that 
ERSI requires to deteet oeeanie features sueh as slieks or internal waves is about 2.5m/s, 
and the same eondition applies to RADARSAT (Staples, personal eommunieation, 2007). 
Both of these results are rules of thumb based on observations. 
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V. RESULTS AND OBSERVATIONS (2) 


This section will test the two methods presented earlier, which retrieve the width 
and amplitude of non linear internal waves using satellite data. Some statistics were made 
to compare the strength of the backscattered signals for the different environments. All 
the following results and observations were collected during the timeframe of ASIAEX 
2001. The different temperature profdes were collected at moorings S4, S5 and S7; the 
satellite images were observed by the Canadian system RADARS AT; and the wind data 
were measured by the scatterometer SSMI when available, or given by the global model 
GFS since no relevant in-situ data were available for this period of time. Within this 
timeframe, only three pictures containing solitary internal waves were recorded for this 
area of interest, and they were all provided by RADARSAT. This system was the only 
SAR system available, as systems such as ERS and JERS experienced some failures. All 
the data were filtered out using an enhanced lee filter available in ENVI to make them a 
little smoother and limit the speckle effects inherent to SAR images. 

Before looking closely at the results obtained, this section will reflect on the 
limitations existing in this analysis. First, the small number of pictures available makes it 
of limited interest regarding statistics, even if it still gives good qualitative information 
about the relationships between solitons, wind and SAR systems. Additionally, all the 
wind measurements used in this analysis were provided by SSMI or GFS, whose 
resolution are respectively 25km and about 100km (1°). The accuracy of this information 
is expected to be around 20° regarding the direction and around 1-2 m/s regarding the 
speed; thus, wind effects must be considered with caution. Finally, one of the most 
difficult parameters to estimate was the speed of the soliton because of the poor vertical 
and horizontal resolutions of the oceanographic array with regard to this project. This 
parameter was, however, essential for the calculation of the width of the observed 
solitons. The speed used to estimate the width of the soliton was a linearly decreasing 
velocity with the bottom, its average value coming directly from the observations at S7, 
S5 and S4 and the satellite images; they were also compared to data provided by Ramp 
(2004). 
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Finally, the main information regarding the characteristics of the soliton or the 
environmental data related to the three SAR images are gathered in Table 3 below. In this 
table, the close and far areas are 5x5 km squares for the smallest and greatest incidence 
angle (no soliton). The soliton area is a 5 x 5km area encompassing the whole soliton 
signature (dark and bright). The dark and bright areas are 1 x 1km squares associated 
with the areas of minimum and maximum caused by the soliton. The mean value 
represents the average brightness within the box. The SAR data were all scaled to bytes 
in that table (0-255). 

Table 3. Summary of the statistics observed for different scenes in the three SAR pictures. 

Day 114 - 1012GMT- wind direction (to)/speed: 25577ni/s - looking direction: 77° 

Lat:21.69°N Processed Range: 113.3km Soliton width: 530m 


Long: 117.22°E Incidence Angle: 27.5° Soliton amplitude: 105m 



Heading: 347° 

Soliton propagation: 300° 


Close 

Far 

Soliton 

Bright 

dark 

Area 

5km X 5km 

5x5 

5x5 

1 X 1 

1 X 1 

Pixels 

160000 

160000 

160000 

6400 

6400 

Mean value 

139.5 

72.2 

90.0 

123.0 

54.2 

Standard 

deviation 

60.0 

44.0 

58.2 

61.8 

41.3 


Day 124 - 1020GMT- wind direction (to)/speed: 210°/5ni/s- looking direction: 77° 

Lat:21.68°N Processed Range: 112.2km Soliton width: 370m 

Long: 117.00°E Incidence Angle: 39.1° Soliton amplitude: 60m 



Heading: 347° 

Soliton propagation: 310° 


Close 

Far 

Soliton 

Bright 

dark 

Area 

5km X 5km 

5x5 

5x5 

1 X 1 

1 X 1 

Pixels 

160000 

160000 

160000 

6400 

6400 

Mean value 

153.4 

107.3 

106.4 

121.4 

87.9 

Standard 

deviation 

62.2 

53.9 

57.4 

60.2 

53.4 


Day 125 - 2215GMT- wind direction (to)/speed: 40°/3ni/s- looking direction: 283° 

Eat:2L87°N Processed Range: 113.6km Soliton width: 470m 

Eong: 116.92°E Incidence Angle: 23.2° Soliton amplitude: 70m 



Heading: 193° 

Soliton propagation: 295° 


Close 

Far 

Soliton 

Bright 

dark 

Area 

5km X 5km 

5x5 

5x5 

1 X 1 

1 X 1 

Pixels 

160000 

160000 

160000 

6400 

6400 

Mean value 

197.7 

72.7 

174.8 

210.8 

147.1 

Standard 

deviation 

49.8 

29.2 

54.1 

51.8 

49.7 
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A. 


FIRST CASE STUDY: APRIL 24, 2001 


Stellite inform Jion 

Heading: 347” 
InddenceAngle: 27.5 
Ascending 



Sale: 
■ 4 -^ 5 Km 


Figure 40. SAR image from RADARSAT on April 24**', 2001 around lOlOGMT as a 
soliton is passing S7. The satellite track, the scale, the observed current 
(moorings) and wind (GFS - l°xl° resolution) are superposed over the picture. A 
zoomed window (3x3 km) is superimposed in the left bottom corner. 



This image shows the presence of a large soliton in the vicinity of mooring S7. Its 
influence is clearly depicted in Figure 41, which represents the temperature observed at 
16 different levels as the satellite passed over the area of interest. Based on these last 
measurements and an approximation of the velocity of the soliton, it is possible to deduce 
the width of the wave and its amplitude using equation 1.7. Further, it was estimated that 
the characteristic half-width was about 550 meters, and the amplitude was observed to be 
95 meters. It must be noticed here that the amplitude may be quite complex to evaluate. 
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Indeed, the amplitude corresponds to the maximum vertical displacement of water as the 
soliton comes through. It obviously varies with the strength of the solution as well as the 
strength of the stratification, which is, in our two-layer model, related to the depth of the 
upper layer. Here the amplitude is assumed to be the relative maximum vertical 
displacement as described by Ramp (2004), and the upper layer was defined by looking 
at the actual density profile. 


SSI vt tlin* - 200 1/1 1 4 



Figure 41. Temperature measurements at S7 before and after the satellite passage. 


To test the methods described by Zheng et ah, different transect profiles were 
observed within the zoomed window in Figure 40. Figure 42 also displays the kind of 
signature obtainable from a SAR image. Such a noisy plot can be confusing and needs to 
be analyzed concomitantly with the picture itself. ENVI offers the opportunity to move 
along the profile and to follow the cursor on the image at the same time. It is then 
possible to identify the most characteristic points in the bright and dark zones related to 
the soliton, and identify the maximum and minimum necessary to apply both methods. 
For this first case, it also appears that the peak-to-peak method presented earlier gives 
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quite good results, sinee the observed eharaeteristie half-width was 514 meters, 
representing a 6 % error. Using the eurve fitting method, the characteristic width is also 
estimated to be between 500 meters and 600 meters as shown in Figure 43. This latter 
figure emphasizes one of the weaknesses of the method, that being that it is difficult to fit 
both the bright and the dark parts of the soliton signature. This problem is inherent to 
SAR systems, since they have both maximum and minimum detection thresholds. From 
this example, the peak-to peak method seems simpler to use and tends to give better 
results. In fact, maximum and minimum can be identified quite easily with an uncertainty 
of one to two pixels, which means less than 50 meters; this resolution is fine enough for 
most of the applications related to solitons. 


Normiiliztd und SAR vulii* Alon9 » tritnjtcl 



Figure 42. Normalized and demeaned transect profile of the SAR information across 
the soliton and in the vicinity of S7. The red crosses represent the significant 
points chosen to determine the characteristic width of the wave, = 678m and 

A = 514m. 
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It must noticed on the SAR image that here the maximum, minimum, bright and 
dark areas are all well-defined, whieh makes the retrieval of the half-width quite easy. In 
the other cases, however, the boundaries between these different areas were not so 
obvious beeause of the wind or the strength of the soliton. 

Best fit of the nornatlizod SAR 5otiton »ind its *th«o r« ti c ail r«pr«5erfetaiti on 



Figure 43. Best fitted eurve (5* order polynomials) and theoretical secant curves 
assoeiated with the observed soliton for different eharaeteristie widths. 


Regarding the amplitude, the results are not as good and the two-layer model 

presented previously barely resolved this parameter. As exposed in equation 1.7, the 

eharaeteristie half-width is greatly dependent on the thiekness of the upper layer, but 

beeause the water eolumn is eontinuously stratified, it is very difficult to set a clear limit 

between the two layers (Figure 44). If one ehooses to set the upper layer to be defined by 

the mixed layer depth or the original depth of the more displaeed water, it appears that 

Zheng’s model tends to underestimate the amplitude of the wave. For example, for a 90 

meter upper layer at S7, the amplitude is expeeted to be only around 30 meters, less than 

a half of the observed value. This part of the model is elearly a limitation to the 

applieation of this theory, and future implementation requires improvements by using 

some stratified models sueh as the one described by Lamb and Yan (1996). Moreover, the 
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secant function in the model makes the amplitude very sensitive to the depth of 
the upper layer, as illustrated in Figure 45. 


Density Profile ft.t 57 



Figure 44. Average density profile at mooring S7 during ASIAEX2001. 



hi 


Figure 45. Sensitivity of the amplitude and the propagation speed of the soliton to the 
definition of the upper layer depth, }\ at mooring S5. The characteristic half¬ 
width is set to be 500 meters. 
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B. 


SECOND CASE STUDY: MAY 4, 2001 


The first case offered the clear depiction of a soliton on a SAR image. However, 
for this second case, the presence of non linear waves was not so obvious, as illustrated in 
Figure 46. These internal waves would have been difficult to obtain without the 
collection of oceanographic data at moorings S5 and S4, because of speckling effects 
caused by SAR. Moreover, even if many solitons were observable on May 4, 2001 in the 
South China Sea, they were all quite small. In this example, only the leading soliton of 
the series will be described. The characteristics of this wave are an upper layer depth 
equal to 80 meters, an amplitude equal to 40 meters and a width around 530 meters . The 
analysis of the current feature utilized an approach identical to that of the first case study, 
and is presented in the following figures. 


Satellite information 

Heading: 347* 
Incidence Angle: 39.1' 
Ascen ding 



Scde: 


5 km 



Figure 46. SAR image from RADARSAT on May 4, 2001 around 1020GMT as a 
soliton is passing S5. The satellite track, the scale, the observed current and wind 
(GFS) are superimposed and a zoom of the area of interest is added in the bottom 

left comer. 
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The presence of solitons in the South China Sea is once again quite obvious in 
Figure 47, which represents the evolution of the temperature profiles for different levels 
at mooring S5. This series of waves is characteristic of a perfect dnoidal group of 
solitons, as presented in the introduction of this paper (Figure 3). Equation 4.9 expresses 
the backscattered energy and shows that it is highly dependent on the depth of the upper 
layer, the width of the soliton and the amplitude of the wave, as well as the wind 
characteristic or the incident angle. For this case, all the parameters which could make the 
signature of the non linear waves noticeable acted against it. Indeed, the incidence angle 
was quite large and greater than 35°, which reduces the backscattered energy. 
Additionally, both the looking direction and the wind direction made the SAR system less 
efficient (Table 3 and 4), the three parameters which characterize the wave tended to 
decrease the backscattered energy. Surprisingly, the signature was still noticeable, and 
these observations may suggest than even weak solitons are observable by SAR images 
as soon as enough wind blows to create capillary waves to return the radar signal. 

SSI ws. tim* - 2001/134 



Figure 47. Temperature measurements at S5 before and after the satellite passage at 8 

different levels. 


75 






















Applying Zheng’s method to this picture was, however, very difficult for this 
case, and the results obtained are not very conclusive. As illustrated in Figures 48 and 49, 
the SAR data observed were quite noisy, and the fit curving shows its limitations, since it 
was impossible to get a clear representation of the presence of the soliton within the 
noise. Thus, the best fit polynomial (fifth order) used here does not at all match the 
theoretical representation of the wave for different widths (Figure 49). Moreover, as 
shown in the zoomed window of Figure 46 compared to Figure 48, the distinction 
between dark and bright areas is not as clear, rendering determination of the location of 
the maximum and minimum more difficult. Indeed, using the peak-to peak method, the 
characteristic width is estimated to be a little bit less than 400 meters, while the 
observations give a result greater than 500 meters. Once again, here the results for the 
amplitude were not very good. 
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Figure 48. Normalized and demeaned transect profile of the SAR information across 
the soliton and in the vicinity of S5. The red crosses represent the significant 
points chosen to determine the characteristic width of the wave, = 520m and 


A = 393m. 
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B«$t fit of fh* normkllitd SKtt solifon knd iff fhtertficail r«o rtftnlsifion 



Figure 49. Best fitted curve (5* order polynomials) and theoretical secant curves 
associated with the observed soliton for different characteristic widths. 


C. THIRD CASE STUDY: MAY 5, 2001 

This last case study offers interesting complimentary information regarding SAR 
detection of non linear internal waves. For this particular picture, the incidence angle at 
the closest edge to the satellite was very small, around 20°, which is a value for which the 
Bragg scattering theory applies less confidently. Secondly, in this case, winds were 
thought to be very light, blowing at speeds close to 3 m/s. This latter value was given by 
the large grid model GFS, and can be considered as doubtful or inaccurate. Nonetheless, 
the value still offers a good order of magnitude, and one can assume that the wind speed 
was within 1 m/s of this value. As shown in Figure 50, the soliton could, however, be 
easily identified for many reasons. First, the orientation of the looking direction with 
respect to the wind direction was favorable. Furthermore, the small incidence angle 
supports a stronger backscattered signal, which, while enhancing the presence of the 
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wave, does not help in diseriminating maximum and minimum. Finally, the 
eharacteristics of the soliton render it strong enough to be detectable by the satellite 
system. The upper layer was estimated to be around 80 meters, while the amplitude of the 
wave at this location was about 50 meters. The width of the soliton was approximately 
300 meters, making it a lot sharper than the two previous waves but with a greater 
amplitude than the second case study. 



Stellite inforrmtior 
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Figure 50. SAR image from RADARSAT on May 5, 2001 around 2215GMT as a 
soliton is passing S7. The satellite track, the scale, the observed current 
(moorings) and wind (GFS -l°xl “resolution) are superposed over the picture. A 
zoomed (3x3 km) window of the area is added in the bottom left comer) of the 

picture. 
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Comparison of the pictures offers interesting insight because the last two cases 
shared close soliton characteristics. In the second case, the soliton was very difficult to 
identify in the noise, whereas the signature is very clear in this last case. The fourth 
section of this chapter will look closely at these observations, while the immediate 
section will focus on the determination of the width and the amplitude using the two 
methods. Once again, only the leading soliton was considered, even though a second and 
a third following wave was noticeable both on Figure 50 and 51. In this last figure, it 
must be noticed that the signature of the soliton is more obvious than in the previous 
case. 


SSI ttm» - 200WI2S 



Figure 51. Temperature measurements at S5 before and after the satellite passage at 8 

different levels. 


From the peak to peak method, the characteristic half-width was calculated to be 
around 450 meters, while the curve fitting method determined it to be between 500 and 
600 meters. These two results are both quite far from the mooring observations, which 
state that the half-width was about 300 meters. Once again, the amplitude was greatly 

underestimated with these two procedures. These complications in the prediction of 

79 
























soliton width most likely arise from the small incidence. The small incidence angle 
creates a strong background reflected energy that probably hid the minimum associated 
with the soliton. The small incidence also tends to make the maximum larger as more 
pixels reach the saturation threshold of the system. As noticeable in the zoomed window 
of Figure 50, it was very difficult to discriminate significant points from the background 
in the expected dark area of the soliton or to set an accurate minimum for the peak to 
peak method. The light wind which blew on May 5 probably did not favor strong 
interactions with the soliton and the appearance of a negative or below normal signal. 
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Figure 52. Normalized and demeaned transect profile of the SAR information across 
the soliton and in the vicinity of S5. The red crosses represent the significant 
points chosen to determine the characteristic width of the wave, = 594m and 

A = 450m. 
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Figure 53. Best fitted curve (5* order polynomials) and theoretical secant curves 
associated with the observed soliton for different characteristic widths. 


D. STATISTICS AND COMMENTS 

In their study, Zheng et ah, promoted their method and claimed that it was 
capable of retrieving the main characteristics of the solitons; this was verified with our 
first case study for the width. The theory behind this elegant model is certainly correct, 
even if some parameters could be described more accurately. However, the last two case 
studies presented above enhance the limitations associated with it. From the three case 
studies, it also appears that the second method based on the distance between minimum 
and maximum is much more conclusive that the one based on curve fitting. The method 
provides an easier way to get the characteristic half-width of a soliton, since it only 
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requires two points. Compared to the eolleetion of huge amounts of oeeanographie data, 
this method ean also be eonsidered aeeurate enough to get the eharaeteristie width 
neeessary for the aeoustie model presented previously or some other applieations 
involving solitons. From these examples, it also seems mueh more eomplieated to get the 
amplitude right, and this part of the theory should be reviewed in a future work. The main 
weaknesses of Zheng’s method were probably related to two important faetors; wind 
eonsiderations and the building of the referenee axes, whieh brings out the eonneetion 
between the wind direetion, the soliton propagation and the looking angle of the satellite. 

1. Wind and Incidence Angle Effects 

As discussed previously, the wind and the incidence angle are factors which can 
significantly affect the detection or characterization of non linear internal waves. The 
following table gathers some information about the three previous SAR pictures and aims 
at comparing the relative effects of the incidence angle, the wind speed and direction and 
the strength of the soliton. Once again, these results must be used carefully, since the 
wind information from the scatterometer SSMI or the model GFS are not very precise. 


Table 4. Relative variation in backscattered energy implied by the presence of a soliton in 


the image along with environmental information. 


Relative 
deviation 
from the 
background 

Bright 

area 

Dark area 

Standard 
deviation 
within the 
soliton 

^(p = (p-(p^l 
cos(A^) 

Wind 

speed 

Soliton 
amplitude 
and width 
(m) 

Day 114 

+46% 

-35% 

+21% 

-178°/-1.0 

7 

105/530 

Day 124 

+7.5% 

-22% 

+2% 

-133°/-0.68 

5 

60/370 

Day 125 

+23% 

-14% 

+9% 

243° / -0.45 

3 

70/470 


The basic results illustrated in this table are the enhanced features brought by the 
presence of a soliton. Thus, for the first case, the brightness is increased by 46% in the 
bright zone and decreased by 35% in the dark, which make the wave very easy to 
discriminate. This result is also observed for the two following cases, though the 
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variations are not as dramatic. In the following two cases, the angle is less favorable 
and the wind is lighter. From only three pictures, it is difficult to determine the most 
efficient wind speed regarding SAR detection of solitons; nor can the importance of wind 
speed relative to wind direction be established. Resolution of these issues requires a more 
thorough study. It seems, however, that winds around 5m/s offer better results than light 
and strong winds. From comparisons of order of magnitudes between the terms of 
equation 4.4, Zheng et ah, stated in their paper that the maximum wind stood by SAR 
systems should be around 5 m/s. Looking at the different cases presented here, one could 
objectively expect this wind threshold to be higher than that and closer to 10 or 15 m/s 
(Raney, 1995). 





2« 24 2» >0 »C 28 «• «2 20 22 2« :• 
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Figure 54. Backscattered energy (and standard deviations) for different incidence 
angles and for three different environments (day 114, day 124 and day 125). 
Statistics associated with the bright and dark signatures of the soliton are 
indicated by blue crosses (1x1 km). The red and green stars indicate the average 
over an area including both the bright and the dark areas. 
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Figures 54 and 55 display the mean backscattered energy and its standard 
deviation for 5km x 5km squares. The discrete values (stars) represent the backscattered 
energy in the bright and dark areas associated with the solitons (1km x 1km) and the 
mean and standard deviation for 5km x 5km areas centered on the solitons. 

These two figures emphasize the effects of the incidence angle on the 
backscattered energy. These figures also clearly show the signature of the soliton. 
Finally, in Figure 55, the comparison between day 114 and day 125 reveals that, for the 
same incidence angle, the returned energy is much stronger in the presence of moderate 
or strong winds. Moreover, it seems that as the angle decreases the system tends to 
saturate much faster. 


Biickf cklttrtd tn«r 9 y obterutd for different incidence )in9le5 And different en oironmentt 
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Figure 55. Backscattered energy observed for different incidence angles for 3 
different environments. The turquoise line represents a scaled scattering 

coefficient ^^^(6*). 
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2 . 


Conclusions and Future Works 


As stated earlier, beeause of the small amount of data available, this study eould 
not aim at giving exhaustive eonelusions about the relationships between SAR systems, 
solitary internal waves and wind. However, this study offered several interesting insights, 
among whieh was the faet that solitons ean be deteeted for low winds around 2 to 3m/s, 
and that the width of a soliton ean be approximated aeeurately enough for many 
operational purposes using Zheng’s method. This study is, however, of a qualitative 
nature, and a thorough experiment should be eondueted to infer more aeeurate 
eonelusions about the wind effeets on SAR deteetion of solitons. The ASIAEX 
experiment was a good start for this kind of analysis, sinee it gathers oeeanographie data, 
aeoustie data and SAR data. The following discussion provides several suggestions for 
improving the current analysis; the expensive cost of such improvements would also be 
debatable, but detailed reflection on this point is outside the scope of this conclusion. 

Clearly, future study requires better spatial resolution. The spacing between S7 
and S5 here was about 25 km and 5 km between S5 and S4, making the observations of 
the propagation speed very difficult. This parameter is, however, essential in the 
determination of the width of the soliton. The timeframe of the whole experiment should 
be increased, much like it was during the Windy Islands Experiment (WISE) in 2005; 
unfortunately, the cost of the SAR images taken during the experiment in question 
precluded their use in this study. The current SAR systems do not repeat their cycles very 
often; to obtain solitons, as well as a large variety of wind conditions and multiple SAR 
images, requires a long experiment, probably around 5 or 6 months. On the other hand, 
studies could proceed with surveillance of many areas at the same time. As stated earlier, 
wind is an important variable in that comparison, and only poor data were available for 
this study. A thorough experiment would require the observations of surface winds by 
buoys concomitantly with the oceanographic measurements. Einally, improvements in the 
method proposed by Zheng et al., are necessary to get correct amplitudes for the solitons, 
as this parameter is even more important than the width, particularly for acoustic 
purposes. 


85 



THIS PAGE INTENTIONALLY LEET BLANK 


86 



VI. CONCLUSION 


Non linear internal waves are eurrently a hot topie in the oeeanographie, 
aeoustieal and military eommunities sinee they eontribute greatly to the water eolumn 
variability. Using data eolleeted during an experiment in South China Sea in 2001, this 
projeet tested the ability of a simple analytie model based on Dyson series to eharaeterize 
aeoustie intensity fluetuations through the retrieval of mode amplitude in the presenee of 
solitons. This theory, presented in the first three ehapters of this report, demonstrates its 
viability by obtaining relevant results for mode 1, one of the most energetie regarding the 
settings of the experiment. Though the results were not as eonelusive for higher modes, 
the model appears promising. The model relies on several assumptions, sueh as the weak 
mode eoupling approximation, the single seattering theory or the simplistie depletion of 
the interaetions between the non linear internal waves and the bottom. Nevertheless, the 
model returns solid qualitative results regarding the mode energy variations. In partieular, 
one of the most satisfying results of this analysis is the predietion of the fast oseillations 
of the aeoustie mode energy as the ISWs move away from the souree. Regarding the 
spreading of the signal observed in the data from ASIAEX 2001, the modeled data did 
not agree with the observations. The modeled spreading was too small, and the inerease 
of the spreading with time predieted by the model was not observable in the in-situ data. 
Other environmental variability, sueh as random linear waves, diurnal effeets or 
mesoseale phenomena, are most likely responsible for the diserepaneies in the model. 

Moreover, it appears that the measurements eolleeted during ASIAEX 2001 did 
not perfeetly mateh those required to eonduet this eomparison. Despite the vast amount of 
data eolleeted during this experiment, it was very diffieult to obtain some of the main 
eharaeteristies of the ISWs, sueh as their amplitude, width and speed, this latter variable 
being a very sensitive input for the mode. In faet, only two moorings were truly available 
to simulate the speed of the soliton and to provide information about the amplitude and 
width. This eonstraint negatively affeeted the aeeuraey of the eomparison. Beeause of this 
drawbaek, it might prove benefieial to eonduet a more thorough analysis of this model 
using experiments with simpler bathymetry, eontinuous measurements and finer spatial 

87 



resolutions, particularly regarding the spacing between moorings and between the source 
and the receiver. Regarding the limitations of this model, future studies should probably 
focus on determining its domain of applicability with respect to water depth and 
frequency. The current study proved that the results calculated by the model were quite 
poor for very shallow waters less than 150 meters, generally due to the complex 
interactions of the waves with the bottom. This model would also most likely break down 
for frequencies greater than lOOOHz, for which the single scattering theory will no longer 
apply. 

As stated in the second and third chapters of this report, the amplitude, 
characteristic half-width and the propagation speed of the soliton proved to be essential 
inputs for the model. This project also described how SAR systems could obtain this 
information. Indeed, over the recent decades, these systems proved their capacity at 
detecting these particular oceanic features and will soon be able to describe their inner 
characteristics. Obviously, the current techniques, such as the one presented by Zheng et 
al., offer only approximate results; however, a few computational and theoretical 
improvements in their methods could significantly improve these methods. This study 
also particularly noted the importance of wind effects in the analysis of SAR images and 
the retrieval of soliton information. This analysis emphasized the evidence of upper and 
lower thresholds regarding wind speed, but the limited amount of data prevents an 
accurate definition. Estimations thus stem from existing data collected by SAR systems 
for winds between 3 and 12 m/s. 
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